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INTRODUCTION 


This report describes the digital computer program SWITCH, which has 
been developed for optimization of multiple-burn vacuum rocket flights* 
The analysis upon which this algorithm is based is described m refer- 
ence [1] and only the equations actually implemented will be described 
here* The report will consist of a description of each subroutine 
used m the algorithm as well as a users guide. 

The SWITCH algorithm is an application of a shooting method to a two- 
pomt boundary value problem of trajectory optimization. The method 
involves guessing initial costate and switching times (separating burn 
arcs from coasting arcs) and propagating state and costate and their 
partial derivatives with respect to initial costate and the switching 
times over a complete trajectory. The partial derivatives are then 
used to * obtain a Newton Raphson adjustment of the initial costate and 
switching times. 

RESTRICTIONS 

The program generates trajectories consisting of alternating coast and 
burn arcs. The first arc may be either a coast or a burn but the last 
arc must always be a burn. All arcs of the trajectories are assumed to 
be positive in time with the exception of an initial coast arc which is 
allowed to be positive or negative in duration (allows starting the first 
burn at any position along initial orbit) . The number of burn and coast 


(1] K. R. Brown, E. F. Harrold, and G. W. Johnson. "Rapid Optimization 
of Multiple^Burn Rocket Flights," NASA Contractors Report CR-1430, 
September, 1969. 



arcs may vary with a maximum of six separate arcs. Discrete changes xn 
vehicle mass are allowed only at the start of coasting arcs , and all 
coasting arcs are assumed to have negative orbital energy (elliptical) . 


The algorithm consists of five basic parts: the input and output rou- 

tines (MAIN and OUT) , the supervisory routine (SWITCH) , the burn routines 
(RKG031, RKSTEP, and YDDRHS ) , the coast routines (COAST and AMULT) and 
the boundary value and linear algebra routines (ADJUST, BVEVA L, and SOLVE) , 
SUBPROGRAMS 

1. MAIN (Input Routine) 

This routine is presented as a sample input routine. All of the data 
necessary to define a case is collected m this routine. The follow- 
ing table defines the input variables m order with FORTRAN read for- 
mats: 

TABLE I - INPUT VARIABLES 


CARD 

VARIABLES 

FORMAT 

1 

NCASE 

12 

2 



3 

PLANS (54) 

18A4 

4 



5 

uk, AMASS, HMAX 

3F10.1 

6 

XOE(I) (I = 1,6) 

6F10.1 

7 

QOE(I) (I = 1,6) 

6F10.1 

8 

C(I) (I = 1,7) 

7F10.1 

9 

LEGMAX, ITMAX 

212 

10 

ATP (1,1); ATP (1,2); 
ATP (1,3); ATP (1,4) 

4F10.1 

11 

ATP (2,1); ATP (2,2); 
ATP (2,3); ATP (2,4) 

4F10.1 


ATP(L,1); ATP(L,2); 
. . .ATP (L,4) 

(~L = LEGMAX + l') 


10 + L 


4F10.1 



CARDS 2 - 10 + L are repeated for each data case. 


The variables are defined on the following page. The mam routine 
also prints the input quantities and initiates a case by calling 
the SWITCH subroutine. 
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VARIABLE NAME 


SYMBOL 


VARIABLE DESCRIPTION 


1) NCASE 

2) PLANS (I) 


3) 

uk 

U 

4) 

AMASS 

m 

o 

5) 

HMAX 


6) 

XOE { I ) 

o 

7) 

QOE{I) 

0 

8) 

LEGMAX 


9) 

ITMAX 


10) 

ATP (1,1) 

T 

1 

ID 

ATP (LEGMAX +1,1) 

T__ 


Number of cases 

ALPHABETIC ARRAY IDENTIFYING EACH 
DATA CASE 

Universal gravitational constant 
Initial mass of the vehicle 

Maximum integration step size 
Initial state vector 

Estimated initial costate vector 

Number of burn and coast arcs 

Maximum number of iterations 
allowed per case 

TIME AT START OF ARC I 
FINAL TIME 


12) ATP (1,2) 

13) ATP (1,3) 

14) ATP (1,4) 

15) C(I) 


Am MASS LOSS AT START OF ARC I. ONLY 

ALLOWED TO OCCUR AFTER BURN ARCS 
(I - 1, LEGMAX) 

m MASS RATE ALONG ARC I 

(I - 1, LEGMAX) 


T 


THRUST ALONG ARC I. IF ZERO ARC 
ASSUMED TO BE COAST (I = 1, LEGMAX) 

Vector of constraints. Definition of 
elements of vector depend on boundary 
value problem under consideration. For 
the rendezvous problem the first six 
components of C contain the state of 
the target at some time which is con- 
tained m C(7). For the 5-constramt 
mission the first three components of 
C contain_ the_ desired angular momentum 
vector (H - R x V) and the next two 
elements contain the first two elements 
of the desired eccentricity vector (E) . 
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2. Subroutine Out (Output Routine) 

This routine controls the output operations for the total program. 

An integer flag variable, NPATH, controls which of three paths the 
program will follow. 

If NPATH is negative, the SWITCH routine has just completed a burn 
arc and the OUT routine prints variables describing the immediately 
preceding coast (if one existed) and burn arcs. The coast arc is 
described by the state and costate at the end of the arc (X and Q) , 
the integral of 1/RADIUS along the arc (PSY) , twice the orbital energy 
(ALPHA) and the calculated coast time (FT) . The burn arc is described 
by the state and costate at the end of the arc (XF and QF) , the mass 
at the end of the arc (AM) and the number of integration steps used 
along the burn arc (ITBURN) . Also outputted at this time are the orbi- 
tal characteristics of the orbit defined by the state of the vehicle 
at the start of the burn arc. The semimajor axis (AAXIS) , the minimum 
radius of the orbit (RMIN) , the maximum radius (RMAX) , the orbital 
energy (ENERGY), the period of the orbit (PERIOD), the angular velocity 
vector (H) and its magnitude (HMAG) , the eccentricity vector (SI) and 
its magnitude (EMAG) , and the radius (RMAG) are used to describe the 
orbit. 

If NPATH is zero, a full iteration has been completed. The resultant 
orbit is described by its orbital characteristics. The total burn time, 
the lengths of the separate burn and coast arcs, the "miss" vector (DC), 
the determinant of the matrix E (partials of constraint functions with 
respect to the initial costate and switching times) , the diagonal elements 
of E, the magnitude of the desired changes m control (DU) and control 
rate scaled by time (DUDT) , the bound CK on adjustments to be made in QO 
and the switching times as well as the vector from whence CK was chosen, 
the requested changes in the initial costate and the switching times, and 
the switching times and initial costate to be used m the next iteration 
are printed. 

If NPATH is positive, the program has completed the case and summary 
tables are printed. The tables include the total burn time, the lengths 
of the separate arcs, the "miss" vector (DC), DU, DUDT, CK, and initial 
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costate for each iteration performed. 


This subprogram is called by SWITCH and calls no other programs. 


3. Subroutine SWITCH 

The SWITCH subroutine supervises the generation of trial trajectories 
and the Newton Raphson iterations. The program accepts as input the 
initial state (XOE) , costate (QOE) , the vehicle description, and 
estimated lengths of the burn and coast arcs. The routine performs 
the necessary initializations and calls the COAST routine to compute 

✓ 

the results of coast arcs and the RKG031 routine to integrate the 

burn arcs. At the end of each arc, a new column of the matrix A 
9 x a 

of partials (*r f — ) is initialized corresponding to the 

OCT ,t , . . . , t 
O 1 1 

time, t., at the end of the arc. 

x 





The parti als are 


f 9x,q 

l 3t 




- (0, 


+ 



m 



0, 0) 


T 

t 

1 


where the sign is negative after coasts and positive after burns. After 
the last arc, always assumed to be a burn, the necessary column of par- 
tials is 

<|f^)| t = <R, R, u, U) T 

F F 

where t is final time. 

F 


The SWITCH routine also calculates the switching conditions and the necessary 
partial s of these conditions with respect to initial costate and the switch- 
ing times. 


The switching condition across a burn arc from t^ to t^ + ^ is 


vw - vv 


- 0 


where 


T v = R T u + iSji 

v IrI 3 


The gradient of the variable, T^, with respect to state and costate is 


3Ty _ ( PU T 3yu T R T .T uR T 

9x 'i' Vi 3 w 5 m 3 


and the partial of the burn switching condition, C , becomes 

3 


3Cr 


3t 


V 


8x,q 


3T, 


V | 3x,q 


3x,q ' \ + l 9 V V ’ ’ * Vl 9X ' q VWz 
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The switching condition across a coast from t^ to t ^ is 


where 


u 


- u 


« 0 


~i41 


|M = (0, 0, -|£p o) 


and the partial of the coast switching condition becomes 


3Cr 

3<I o' fc l"’ 



3 juj | 3x,q 

3x,q 't 1+1 3q o( t l( 


. t 


i+l 


3 jul | 3x x q 

9x,q t 3q , t_ , , . . ,t 
^ i ^o 1 1 


After the trajectory has been calculated SWITCH calls BVEVAL to evaluate 
the six terminal constraints and then calls ADJUST to calculate the Newton 
iterate desired. Another iteration is then performed with the new costate, Q01, 
and switching times (ATP (1,1)). 


At the end of each iteration a test is made to see if the maximum number 
of iterations (KMAX) has been exceeded or if the new adjustments are less 
than the tolerances which indicates convergence. When either of these occur 
NO is set equal to zero and a final trajectory is computed propagating 
only the state and costate (not the parti als) . 


The SWITCH program also calls OUT (the output routine) and controls the 
amount of printout. 
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4 * Subroutine RKG 031 


This program propagates state , costate , and their transition matrix 
across a burn arc, using a numerical integration method especially 
adapted for the purpose. The inputs accepted as explicit arguments 
are the initial state and costate XO and QO and the initial matrix of 
partials Z, together with a desired error level indicator EVT, a step 

A 

size limit HMAX , the number LEGMAX of the final leg, and the indicator 
NO which, if zero, suppresses computation of the partial derivatives* 

Other inputs accepted via COMMON are the number LEG of the current leg, 
and mass, thrust, and timing for that leg m the ATP array. Explicit 
outputs are state and costate XF and QF and the matrix of partials Z at 
the end of the arc* Also an integration step count ITBURN is an output 
through COMMON. 

The program operates by repeatedly calling on RKSTEP to perform Runge- 
Kutta numerical integration steps* For the sake of efficiency, these 
steps are grouped into combination steps consisting of three short steps 
of integration of state and costate equations only and one overlapping 
triple length step of integration of state and costate equations and 
variational equations* Thus the computationally expensive right hand 
sides of the variational equations are called only 25% as often as the 
state and costate equations. This is permissible since the matrix of 
partials is relatively non-critical for accuracy, for it is used only m 
linearized Newton-Raphson adjustments which are only approximate anyway* 
Also the redundancy m integration of state and costate is useful m esti- 
mating integration accuracy so that step sizes can be continually reopti- 
mized* The error estimate is obtained by the Richardson extrapolation 
technique, which (m this case) estimates the error m the result of the 





three short steps to be one eightieth the difference between the 
result of the three short steps and the single long step. This 
error estimate is added m to eliminate the error to the extent 
that it was correctly estimated, and it is also used to determine 
the best step size for the next combination step. The desired 
error level for one step is assumed to be equal to the length of 
the step divided by the length of the entire burn arc, times a 
desired level for the total error over the burn arc. Estimated 
error is compared with desired error after each step, and step 
size for the next step is then altered by the fourth root of the 
ratio between the two. [Note: actual adjustments involve error 

squared so that an eighth root actually appears m the program.] 

When the burn arc is the last of the trajectory, RKG031 also calls 
YDDRHS to obtain partials of state and costate with respect to final 
time TF and stores these partials as a final column of the Z matrix. 

5. Subroutine »RKSTEP 

This subroutine performs a single fourth-order Runge-Kutta integration 
step on a system of second order differential equations for a matrix 
function Y consisting of 6 rows and at most 13 columns. The calling 
program supplies the following input values as explicit arguments: 

YN = initial value of Y 

YDN - initial value of the first derivative of Y 
TN = initial time 
H = step size 
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N = flag variable 


LEG = number of columns in excess of 6 to be advanced as 
explicit arguments 

and RKSTEP calculates the following values which are returned: 

YNl = value of Y at TN + H 
YDNl = value of the first derivative of Y at TN + H 
If the flag variable, N, is positive, then the first 6 + LEG columns 
of Y are advanced. Otherwise, only the first column is advanced. 

The method of attack consists of letting Y = F (Y,t) represent the 
differential equation of the matrix Y. Then, with 
= H*F(YN,TN) 

D 2 - H*F (YN + |(YDN + ) , TN + |) 

D 3 = H*F (YN + H* (YDN + ) , TN + H) 

YN1 and YDNl may be evaluated as 

YNl = YN + H* (YDN + -^(D, + 2*D„) ) 

6 1 2 

YDNl = YDN = !-(D-, + 4*D„ + D) 
o 1 2 o 

Evaluations of the form H*F(Y,t) which appears in the above equations 
are accomplished by calls to subroutine YDDRHS. 
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6 * * Subroutine YDDRHS 


This subroutine performs an explicit evaluation of a matrix function 
YDD = H*F(Y,T) consisting of 6 rows and at most 13 columns. The calling 
program supplies values of the input variables Y,H and T and a flag vari- 
able N as explicit arguments. YDDRHS then calculates an appropriate value of 
YDD which is returned as an explicit argument. If the flag variable, N r is 
positive, then the first 6 + LEG columns of YDD are computed. If N is non- 
positive, then only the first column of YDD is calculated. The parameters 
tfK, AM, LEG and ITBURN and the matrix ATP are transmitted via labeled common 
block /CCPINJ/. 


The matrix YDD may be partitioned 


YDD « 


- wdd 

1 UDD 


and Y may be partitioned 


I R * T7 

t 5 1 J 


where HDD, UDD, R & U are 3x1 matrices (3-vectors) and WDD and W are 
6 x (LEG +5) matrices. Using standard vector-matrix notation, RDD, UDD, 
and TOD may be expressed m terms of R, U and W as 


RDD = - 


UR 


ATP (LEG ,4) U 

3 AMI ]uT 


UDD = E _ J4L 


WDD = BW 
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where }i = UK and the matrix B is given by 


B 





— 


B_ , 

l 

B 


_ _11 

\ 

- - 1 ? - 


B 

~ 3 

B 


21 

t 

11 



• 



The 3x3 sub-matrices of which B is composed are given m terms of 
R and U as 


b„ - -JUs-Sy . ! . 

11 Ir ! 3 IrI 2 


B 12 " (I 3 - 


UU 1 ATP (LEG/ 4) 
|UJ2 J M AMI 


T T 

r RR , 3yR ,_„T . „_T, 3y 

B 21 = (I 3 ” 5 e" + RU + UR ) - 


l R l ! r 
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7 . Subroutine COAST 


This program takes the initial state (RO,VO) , initial costate 

(DRO,DVO) , plus the desired coasting time and calculates the 

final state (R,V) , final costate (DR,DV) , and the matrix of par- 

* 

tials of final state and costate with respect to initial state 
and costate. The state and costate vectors are transferred by 
labeled common ACOAST. The calling sequence contains as inputs 
the gravitational constant UK, a flag variable NO, ATP and LEG 
which are used to calculate the coasting time. IF NO equals zero 
only state and cost ate are propagated. IF NO is nonzero the par- 
tials are also calculated. The calling sequence also includes PHI 
and DP HI which are (6x6) matrices of partials. PHI is the par- 
tials of final state (costate) with respect to initial state (co- 
state) . DPHI is the partials of final costate with respect to 
initial state. The labeled common, WCOAST, transfers PSY (the in- 
tegral of |r| ^ over the coast) , ALPHA (twice the orbital energy), 
and FT (the calculated coast time) to the output routine (OUT) . 

The general solutions for the two-body problem are presented m ref- 
erence [2] . The equations presented in the paper by Goodyear are im- 
plemented as presented for the closed form calculation of R, V and PHI. 

It was noticed that the partials of final costate with respect to 
initial state (DPHI) was equal to the formal differential of the PHI 
matrix corresponding to letting costate be considered to be the dif- 
ferential of state (this fact is derived [1] ) . The variable names m- 

[2] W. H. Goodyear. ''Completely General Closed-Form Solution for Coor- 
dinates and Partial Derivatives of the Two-Body Problem, 1 ' The Astro- 
nomical Journal , Vol. 70, No. 3, April, 1965, pp. 189-192. 
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ternal to COAST are a one-to-one mapping of the names used m the 
Goodyear paper. The differential of a variable is specified by a 
D in front of the variable name (,e.g. SO -»■ DSO) . It should be 
noted that final costate is evaluated as the differential of final 
state (DR,DV) , 

This program is called by subroutine SWITCH and calls a special 
purpose matrix multiply routine AMULT, 
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VARIABLE DEFINITIONS 


VARIABLE NAME 

COAST ROUTINE VARIABLE DESCRIPTION 


1) 

RO 

Initial radius 

2) 

VO 

Initial velocity 

3) 

DRO 

Initial control (u) 

4) 

DVO 

Initial control rate (u) 

5) 

R 

FINAL RADIUS 

6) 

V 

FINAL velocity 

7) 

DR 

FINAL CONTROL 

8) 

DV 

FINAL CONTROL RATE 

9) 

PHI 

PARTIALS OF FINAL STATE/costate 
with respect to initial STATE/costate 

10) 

DPHI 

PARTIALS OF FINAL COSTATE 
with respect to initial state 

ID 

PSY 

as defined in reference [^L] (Goodyear) 

12) 

ALPHA 

Twice orbital energy 

13) 

FT 

Calculated coast time 
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8 . Subroutine BVEVAL 


This program is used for evaluation of the six right end and -(PER- 
HAPS) transversality conditions which are needed to specify the 
desired final orbit. The program computes the difference DC be- 
tween the desired value of the terminal constraints and their 
actual values as determined from XF and QF. It also computes the 

matrix E of partial derivatives of the right end conditions with ^ 

" ! * 

respect to QO and the switching times. Two versions of subprogram 
BVEVAL are presented here: one for a five-constraint orbital in- 

fection with time origin free, and the other for rendezvous which 
involves a full complement of six constraints on final state. 

The five-constraint version of the program assumes the six terminal 
conditions consist of the angular velocity vector, the first two 
components of the eccentricity vector, and one transversality variable 
at T . The program calculates the partial derivatives of these con- 
straints with respect to XF and QF and multiplies them by the matrix Z 
(partial derivatives of XF and QF with respect to QO and the switching 
times) to form the matrix E of partial derivatives of the terminal 
constraints with respect to QO and the switching times. 

The rendezvous version of the program uses as constrained terminal vari- 
ables the full state vector at final time. In this case the C vector 
which specifies the terminal manifold contains the state of the target 
yehicle in its first six components and the corresponding time in its 
seventh component. In order to form the DC vector, the state of the 
target vehicle is propagated to and subtracted from the actual XF. 

The matrix of partials of right end conditions with respect to 00 and 
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the switching times becomes for this case the first six rows of the 
* 

Z matrix except that the partial derivative of target state with 

respect to T is subtracted from the last column of E to produce 

the matrix of partial derivatives of the DC vector with respect 

to QO and the switching times. The rendezvous version of BVEVAL 

prints the state of the target vehicle at T following each iteration and 

F 

prints the orbital characteristics of the target orbit after the first 
iteration * 

The BVEVAL subroutine is called from SWITCH and calls no other routines* 
9* Subroutine ADJUST 

This subprogram adjusts the initial costate vector (QO) and the 
switching times (ATP (1,1) ) based on linearized relations between 
these quantities and a vector C of constraint functions which is 
desired to change by an amount DC (the "miss 11 vector) . The inputs 
are QO, ATP, E (the matrix of partials of C with respect to QO and 
the switching times) , Z (the matrix of partials of state and costate 
with respect to QO and switching times) and LEGMAX. The vector DC 
appears as an additional column of the E matrix* 

The explicit outputs from the routine are new initial costate (QO) 
and new switching times (ATP (1,1)), the magnitude of the desired 
change in the initial control vector (DU) and the magnitude of the 
desired change m the initial control rate vector times total mission 
time (DUDT) * The bound CK on adjustments allowed on QO and the switch- 
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mg times and the vector A containing the maximum allowable changes 
in the switching times and QO are passed directly to OUT via labeled 
common WADJ. The matrix E may also be considered an output since it 
is altered during the execution of this program. 


The desired change m the 00 vector and the switching times (t^; 
i = 1, LEGMAX) is obtained by solving the following simultaneous 


equations ; 


Ac = [ 


3c 


3QO , t^ , • • • / 1^ 


aqo 

At. 


At 


by calling subroutine SOLVE. The desired change in the QO and the 
t^'s is placed in the LEGMAX + 7 column of the E matrix. 


In order to insure convergence only a fraction ^CK^of the desired 
linearized changes m initial costate and switching times is ac- 
cepted. CK is determined as the largest number (SI) compatible 
with the following restrictions: the magnitude of the change m 

the initial control vector (DU) is limited to .2, the magnitude 
of the change in the initial control rate (DUD) is limited to .0003 
and the length of the burn and coast arcs is allowed to change by 
at most 50%. [This prevents arc durations from disappearing or 
becoming negative.] However, no restriction is placed on the allow- 
able change in the duration of initial coasts. This allows for the 
possibility of negative initial coast durations, i.e. the first burn 
may begin before the (arbitrary) epoch t Q if it is optimal for it to 
do so. 
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The new QQ and switching times are then defined by adding to their 
previous values increments equal to CK times the Newton Raphson 
increments. The new QO vector is then normalized such that its 
first three components (u) have a magnitude of unity. 

This subprogram is called by SWITCH and calls SOLVE. 
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10 . Subroutine SOLVE 


This program solves a system 
6+LEGMAX 

A X = A ,i = 1 , . . .6 + LEGMAX; L = 7 + LEQ4AX 
lj 3 lL 

of linear simultaneous equations for the unknowns (X^ ) . The inputs 
are the elements of the matrix A, with last column (A (I, 7 + LEGMAX)) 
containing the right hand sides of the equations and LEGMAX specify- 
ing the number of equations to be solved. The output is the same 
matrix A with the last column replaced by the solution vector X and 
the other columns destroyed (by the Gauss-Jordan reduction transforma- 
tions used to solve the linear system) . 


I 

1=1 


The solution is obtained directly by Gauss-Jordan reduction without 
generating the inverse of the coefficient matrix as an intermediate 
result. In order to avoid undue round off errors, at the ith stage of 
the reduction process the unreduced rows are sorted such that element 
l + 1 of row i+l is at least as large m magnitude as the corresponding 
element of any later row. This has no effect on the ordering of the solu- 
tion vector. 


This program is called from ADJUST and calls no other routines. 
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APPENDIX I 


SWITCH Program Listing 



FILE:< SWITGHD FORTRAN ?1 


CAMBRIDGE MONITOR SYSTEM 


******** MAIN ROUTINE **** ****** 

IMPLICIT REAL*3 (A-H,0-Z) 

COMHON/CCPIN J/UK , AT? (7,4) , AM, LEG, IT BURN 
COMMON/WIN/PLANS (54) , NOUT, LOGIC 
DIMENSION C (12) , QOE (6) ,X0E (6) 

FORM AT (7F10. 1) 

NCASE , RIGHT ADJUSTED IN COLUMNS 1-2 
CASES FOLLOWING 


C READ INTEGER VARIABLE 
C IS THE NUMBER OF DATA 
NIN=2 
NO UT=8 

READ (NIN, 3) NCASE 
NC=0 

5 N C=N C+ 1 
LOGIC-O 

WRITE (NOUT, 6) NC 

6 FORMAT (1 HI ,13H CASE NUMBER ,12,//) 
READ(NIN, 10) PLANS 

10 FORMAT (18 A4) 

WRITE (NOUT, 10) PLANS 

C FIRST THREE CARDS OF EACH CASS CONTAIN 
C COLUMNS 1-72 DESCRIPING TH3 CASE BEING 
READ (NIN, 1 ) UK , Ail AS S , H d A X 
THE NEXT CARD CONTAINS THE GRAVIATIONAL 
THE INITIAL MASS OF THE VEHICLE , AMASS 
A LL0NA3LE INTEGRATION STEP SIZE , HMAX 


WHICH 


ALPHAB ETIC 
SOLVED 


INFORMATION IN 


CONSTANT , UK 
, (COLS 11-20) 
(COLS 21-30) 


FORMAT FI 0.1 


TFI SI< VECTOR X OE WHICH IS THE INITIAL STATE 
5F10.1 FORMAT. 


READ (NIN, 1) XOE 
THE NEXT CARD CONTAINS 
OF THE VEHICLE UNDER A 
READ (NIN, 1 ) QOE 

NEXT IS THE INITIAL COSTATE VECTOR , CONTROL AND CONTROL RATE (QOE) , 

UNDER 6F10.1 FORMAT. 

READ (NIN, 1) (C (I), 1=1, 7) 

THE 7-VECTOR , C , IS READ NEXT UNDER A 7F10.1 TORMAT. THE VECTOR 
CONTAINS THE DESIRED VALUES OF THE RIGHT END CONSTRAINTS OF THE 
MISSON. THE BOUNDARY CONDITIONS FOR A RENDEZVOUS MISSON ARE THE STATE 
(LOCATED IN FIRST SIX COMPONENTS OF C) AND CORRESPONDING TIME (IN 
DEFINING THE ORBITAL POSITION C? THE TARGET VEHICLE. FOR THE 
5-CCKSTRAINT MISSON THE FIRST THREE COMPONENTS OF C CONTAIN THE 
DESIRED ANGULAR MOMENTUM VECTOR (H= R X V) AND THE NEXT TWO ELEMENTS 
OF C CONTAIN THE FIRST TWO ELEMENTS OF THE ECCENTRICITY VECTOR (E) . 

READ (NIN, 3) LEG MAX, I MAX 
THE NUMBER OF BURN AND COAST ARCS , 

ITERATIONS ALLOWABLE , TMAX /APPEAR 
i FORMAT (212) 

FORMAT (4F10.1) 

Ll=LEGllAX + 1 

READ (NIN ,4) (ATP (1,1) , ATP (1,2) , ATP (1,3) , ATP (1,4) 

MATRIX CONTAINS A ROW DESCRIPING VEHICLE CHARACTERISTICS FOR EACH 
SEPERA7E ARC. ATP (1,1) - TIMS AT START OF ARC I , ATP (1,2) - MASS LOSSS NlOOiiSO 

AT START OF ARC I (ONLY AFTER 3 J BN ARCS) , ATP (1 ,3) - MASS RATE ON ARCSWI00500 
I , ATP (1 ,4) - THRUST ALONG ARC I (IF ZERO ASSUMED COAST). ATP ARRAY SWICQ510 
ALSO HAS A ROW CONTAINING FINAL TIME (TF=ATP (LEGMAX+ 1,1)). THE MATRIX SWI0O520 
IS INPUTTED ONE ROW PER CARD UNDER A 4F10.1 FORMAT. SKI00530 

WRITE (NOUT, 100) UK, AM ASS , HM AX , IMAX , LEG MAX SWI00540 


LEGMAX 
ON THE 


, AND THE 
NEXT CARD 


C ATP 


SWI00010 
SWI 00020 
SWI00030 
SWI00040 
SWI00050 
SWI00060 
SWI00070 
SWI00080 
SWT 0009 0 
SWI00100 
SWI001 10 
SWI00120 
SWI 001 3 0 
SWI00140 
SWI00150 
SWI00160 
SWI 001 70 
SWIQO 180 
SWI00190 
SWI00200 
SWI 00 21 0 
SWI00220 
SWI00230 
SWI002UO 
SWI00250 
SWI0026 0 
SWI 00270 
SWI00280 
SWI00290 
SWI00300 
SWIOO 310 
SWI00320 
SWI00330 
SWI0034 0 
SWIOO 350 
C (7) ) SWIC036U 
SWI 0037 0 
SWIOO 380 
SWI0039Q 
SWI00400 
SWIOO 4 1 0 

MAXIMUM NUMBER OFSWI00420 
WITH A 212 FORMATSWI0043 0 

SWl00ti4Q 


, (COLS 1-10} , 
AND THE MAXIMUM 
ALL UNDER 


[= 1 , LI) 


SWIOO 450 
SWIOO 460 
SWI0047 0 
SWIOO 480 
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00 FORM AT { ’ OGR AVITATIONAL CONSTANT^ 'FI 0 .2 , 1 INITIAL MASS= 
1 INTEGRATION STEP=* F 1 0 . 3 ,/ * MAX NUMBER OF ITERATION S=’ 
20F SFPERATE ARCS=»,I3,/) 

WRITE (NOUT, 10 1) (ATP (I, 1) , ATP (1,2) , ATP (1,3) ,ATP (I, ft) ,1= 

01 FORMAT (» O' 35X 1 ATP ARRAY '// ,3X ’TIME 1 1 3X r DELTA MASS 1 10 X» 
1 13X'THBUST'//(1X,4E20.8} ) 

WRITE (NOUT, 102) XOE,Q0E, (C(I) ,1=1,7) 

102 FORMAT (17H0INITIAL STATE XO , 6? 1 7. 7/1 6H ESTIMATED QO 
116E0DESIRSD FINAL C,7E16.8/) 

CALL SWITCH (A MASS , LSGH A X, I.1AX, C,Q0 E,X0 E, HMAX) 

IF (N CASE . GT . NC) GO TO 5 

STOP 

END 


* El 6. 8 , 1 MAX SWI 0055 0 
13,' NUMBER SWI00560 
SWI00570 
1 , LI) SWI0058Q 

MASS RATE* SWI00590 


,6F17 . 7/ 


SUBROUTINE OUT ( X, Q , XF ,Q F , L EGMAX , NPATH) 

IMPLICIT REAL-8 (A-H,0-Z) 

C ROUTINE FOR PRINTING DATA. X , Q -STATE AND COSTATE AT BEGINNING OF 
C ARC JUST COMPLETED. XF , Q? - STATE AND COSTATE AT END OF ARC. NPATH 
C INTEGER VARIABLE CONTROLLING PRINT OPTIONS { IF NEGATIVE - END OF LEG 
C IF ZERO - END OF ITERATION , IF POSTIVS - END OF CAS 5 ) . 

COMHON/CCPIN J/UK , AT? (7,4) , A M , LEG ,1 TBUPN 
COMMON/ WI N/PL ANS (54) , NOUT, LOGIC 
COMriON/W AD J/A ( 8) ,C< 

COilMON/WCO AST/PSY, ALPHA , FT 

COMM ON/ WS WIT/3 (12,13) ,DC (12) ,D U ,D UD T, S VT , BUR NT , KCOUNT 
DIMENSION X (6) ,H(3) , El (3) ,TIME (50,7) , DC1 (50, 12) , DU1 (50) ,DUDl (50) 
DIMENSION KC 1 (50) , Q (6) ,XF(6) ,QF{6) ,CK1 (50) ,CQ0 (50,6) 

IF (NPATH) 1,2,2 

1 LEG1 =LEG-1 

IF (LEG . LT . 3) WRITE (NOUT ,40) KCOUNT 
40 FORMAT (1 HO , 17HITE RATION NUMBER ,13,/) 

IF (LEG. EQ. 1) GO TO 2 
WRITE (NOUT, 10) LEG1 ,X,Q 
C OUTPUT STATE AND COSTATS AT END OF COAST 
10 FORMAT ( 20X , 9 HCOAST ARC ,/25X ,4HLEG=, 

112, 1X12HSTATE AT END , 1 X , 6E 14. 6 , /30X, 1 4HCOST ATE AT END, 1 X , 6E 1 4 . 6) 
WRITE (NOUT, 22) PSY, ALPHA, FI 

C COAST ARC PARAMETERS - PSY = INTEGRAL OF 1/R OVER COAST , ALPHA = 

C TWICE THE ORBITAL ENERGY , FT = CALCULATED COAST TIME 
22 FORMAT (30X, 4HPSY=, E 1 4 . 6 , 5 X, SHALPH A= ,E 1 4 . 6 , 5X, 22HC ALCUL ATED COAST 

1IME= ,E14. 6) 

2 B (1)=X (2)*X(6) -X (3) (5) 

C CALCULATION OF THE ORBIIAL ELEMENTS FOR THE ORBIT DEFINED BY X 
H (2) =X (3) *X (4) -X (1) *X (6) 

H(3) =X(1) *X (5) -X (2) -X (4} 

HM=DSQRT (H (1) **2+H (2) **2 + H (3) **2) 

RM=DSQRT (X (1 ) **2+X (2) **2+X(3) * + 2) 

El ( 1 > — - (X ( 1) /KM+ (H (2) *X (6) -H (3) *X(5) ) /UK) 

El (2)-=- (X (2) /RM+ (H (3 ) * X (4)-H{1) -X ( 6) ) /UK) 

El (3) =- (X (3) /RM+ (H (1) -X (5) -H (2) *X (4) )/UK) 

EM=DSQRT (El (1) --2+S1 (2) *-2-*-E1 (3) *-2) 

ENESGY=— UK/PfH- . 5* (X (4)**2+X (5) **2+X (6) **2) 

AAXIS=-UK/ (2 .0-ENERGY) 
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RMIS=AAXIS* (1 .O-EM) SWIOI 07 G 

RMAX=AAXIS*(1 .0+EM) SWI01030 

PERIOD= (6. 2831853) *DSQRT { D A 8S ( AAXIS**3/UK) ) SWI01090 

IF (NPATri. GE.O) WRIT'S (NOUT ,41) SWI01100 

41 FORMAT (22X,30HRESULTANT ORBIT SPECIFICATIONS) SWIOI 110 

WRITE (NOUT ,11) AAXTS , 5MIN, EM AX, EN ERGY ,P ER IOD, FM, H,EH,E1 , RH SWI01 120 

11 FORMAT(25X, 15 HSEMIfl A JOE AXI3=,E 1 4 .6 , 1 X5HRMI N= ,E14 . 6 , 1X5HRMAX=, El 4 SWI01130 

1 . 6, 1X7HEN3RGY=, El 4 . 6 ,/25 X, 7 HPERIOD=, E 14 .6, 1 X5HHMAG=,E14 . 6 , 1 X8HH VESWI01 140 
2CT0R, 3E 1 4. 6 ,/25x5HEM AG=, 31 4 . 6 , 1 X8HE VECTOR, 3E1 4. 6 , 1 X5HRMAG= , E14 . 6) SWIOI 150 

IF(NPATri) 5,6, 6 SWI01160 

5 WRITE (NOUT, 12) LEG, X? ,QF ,A M, IT3 U RN SWI01170 

C OUTPUT STATE , COSTATE AND MASS AT END OF BORN SWI01130 

12 FORMAT (20X,8H30RN ARC , /25X , 4HLEG=, 12,1X1 2H3T ATE AT END , 1 X6E1 4 . 6 , /SWIOI 19 0 
13 OX, 14HCOSTATE AT END , 1 X6S 1 4 . 6 , /3 OX , 1 9HI1 ASS AT END OF LEG=,E14.6, SWI01200 
25X, * NUMBER CF INTEGRATION STEPS IN BURN ARC=*I4) 

IF (LEG.LT.LEGMAX) RETURN 
KC=MOD ( KCOUNT- 1 , 50) 

KC=KC+1 
DO 7 1=1, LEG 

7 TIME (SC, I)=ATP {1+1,1) -ATP (1,1) 

RETURN 

6 L6=6+LEG 

TIME (KC , 7) =BU RN T 

WRITE (NOUT ,13} TIME (KC,7) , (TIME(KC,I) ,1=1, LEG) 

C OUTPUT TOTAL BURN TIME AND LENOIR Or ARCS 


13 FORMAT (20X, 16HTOTAL EUSN T I M E=, E 1 4 . 6 , 1X9HARC TIMES , 1X4S14 . 6 , /6 1 X 
12E14.6) 

IF { NP ATH . GT. 0) GO TO 3 
DETE=1 .0 
DO 4 1= 1 , L 6 
u DETE=DETE*S (1,1) 

WRITE (NO UT ,14) DC ,DBTE , (3 (1,1) ,1=1, L6) 

PRINT VECTOR OF ERRORS IN 300 N CARY CONDITIONS (DC) AND DETERMINANT AND 
DIAGONAL ELEMENTS OF MATRIX OF P AKTIALS OF BOUNDARY CONDITIONS WITH 
RESPECT TO INITIAL COST ATE AND SWITCHING TIMES 

14 FORMAT (25X,2HCC,1X6E14. 6 ,/2SX,6El4. 6,/25X , 1 7HDETERMINANT OF S= , 
1E14. 6,1X1 3HDIAGO NAL OF E , 4E 1 4 . 6 ,/25X, 8E 1 4 . 6) 

Li = LEG + 1 
L2=L 1+1 

WRITE (NOUT, 15) DU, D UDI,CK ,EVT , (A (I) ,I=1,L2) 

15 FORMAT (25X,3HDU=,E14. 6, 5H DUET=, E14.6 ,4 H CK=, E14. 6, 5H EVT=,E14.6/ 
125X, 1 0HCK=MIN OF ,7E13.5) 

L7=LEG+7 

WRITE (NOUT, 30) (E(I,L7) ,1=1, L6) 

30 FORMAT (/25X66HCHANG5 REQUESTED IN INITIAL COSTATE , SWITCHING TIMES 
1 AND FINAL TIME,/25X , 6H 1 4 . 6 , /35X , 6E1 4 . 6) 

WRITE (NOUT, 16) KCOUNT, QF, (ATP (I, 1) ,1=2, Li) 

16 FORMAT (IX, 24FEND OF ITERATION NUMBER ,I3,/1 5X,7HNE? QO ,6E14.6,/ 

1 1 5X , 1 6 H N C W SWITCH TIMES , 1X5E14 . 6 , //) 

DO 3 1=1,6 
CQO(KC,I)=QF(T) 

DC 1 (KC,I) =DC (I) 

8 DC1 (KC,I+6)=DC (1+6) 

CK 1 (KC) =CK 
DU1 (SC) =DU 


SKI 01 21 0 
SWIOI 220 
SWI01230 
SWI01240 
SWI01250 
SWI01260 
SWI01270 
SWI 01 28 0 
SWI01290 
SWI 01 300 
SWI 01 31 0 
SWIG 1 320 
SWI01330 
SWI01340 
SWI 01 350 
SWI01 360 
SWI0137 0 
SWI01380 
SWI01 390 
SPI0140G 
SWI01 410 
SWIM 420 
SWI01430 
SWI 01 44 0 
SWI 01 450 
SWI01460 
SWIM 470 
SWI01480 
SWI 01 490 
SWIO 1 500 
SWIC1510 
SWI01520 
SWI 01530 
SFI01540 
SWI0155O 
SWIOI 560 
SWI 01 570 
SWIOI 580 
SWI0159O 
SWI01600 
SWIOI 61 0 
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DUD1 (KC) = DUDT 
KC1 (KC) =KCOUNT 
RETURN 
KCP =K COUNT 
IF (KCP. GT. 50) KCP=5 3 
KC 1 (KC) =KCOUNT 
WRITE (NOUT, 11 5) PLANS 
FORMAT (1 HI , (1 8A4) ) 

WRITE (NOUT, 17) 

FORMAT (1H0 ,50X, m-lSUMMARY TABLES,//, 1X9 HIT ERATIO S, 1 X 1 5HTOT AL 1 
1 TIME, 21 X, 29 HL ENG TH OF 3 URN AND COAST ARCS ,/2 X6 HNUdBSR ,//) 

DO 9 1=1, KCP 

WRITS (NOUT, 18) KC 1 (I) ,TIME (1,7) , (TIME (I,J) ,J=1,LEG) 

FORMAT (3X,I3,5X,7B15.7) 

WRITE (NOUT, 19) 

FORMAT ( 1H0,//1X 9 HITS RATION, 39X,36H ERROR IN BOUNDARY CONDITIONS- 
11-8) ,/2X,6H NUMBER,//) 

KC1 (5C)=KC1 (KC- 1) + 1 
IF (KCOUNT.LE. 50) KCP=KC?-1 

WRITE (NOUT, 20) (KC1 (I) , (DC1 (I,J) , J=1,8) ,1=1, KCP) 

FORM AT (3X13, 5X,8E1 4. 6) 

WRITE (NOUT, 21) 

FORMAT (1 HO, //IX 9HITS RATION ,2oX,8HDC (9~ 12) , 33X,2HDU, 12X, 4HDUDT, ' 
1 HCK ,/2X6 hNU MBER ,//) 

WRITE (NOUT, 18) (KC1 (I) , (DC1 (I,J) ,J=9,12) ,DU1 (I) ,DUD1 (I) ,CK1 (I) ,! 
1 , KCP) 

WRITE (NOUT ,51) 

FORM AT ( 1 HO , // 1X9H ITER ATION, BOX, 37HNEW QO GENERATED 3Y PRESENT ! 
1 ATI ON, /2 X6H NUMBER ,//) 

J=0 

WRITE (NOUT ,52) J, Q 
FORMAT (3X13, 5X, 63 16. 8) 

WRITE (NOUT, 52) (KC1 (I) , (CQ9 (I, J) , J=1, 6) ,1=1 , KCP) 

RETURN 

END 


SWI01959 


SUBROUTINE SWITCH (A’l ASS , LEGM AX, IMAX ,C , QOE ,X0 E ,H MAX) 
IMPLICIT REALMS (A— H,0-Z) 

COMMON/CCPINJ/UK, AT? (7,4) , A M , LEG, ITBU RN 
COMiiON/WSWIT/E ( 12,1 3) ,DC(12) ,DU , D UDT , E VT, BURNT , KCOU NT 
COM MOK/ACOAST/XO (6) , Q9 (5) ,XP (6) , QF(6) 

DIMENSION C (12) , Z ( 1 2 , 1 2) , X 03 (6) , Q01 (6 ) 

DIMENSION QOE (6) ,DUMMY(13) ,?HI{6,6) ,DPHI(6,6) ,DZ(12,12) 
KCOUNT=0 
NO= 1 

L7=7-*-LEGMA X 
L6= 6+LSGMAX 
E V T= 1 . E-8 
KHAX=IMAX 
L1=LEGHAX+1 
DO 2 1=1,6 
Q01 (I) =Q0E (I) 

J. «!. U J.X JLV J.J. ^4. J, 4/ 4* 4. J. vV 4. 4. 4 


SWI00010 
SWI00020 
SWI00030 
SWI00040 
SVJI00050 
SWI00060 
SWIQ0070 
SWIG9080 
SWI00090 
SWI001 00 
SUI00 110 
SWIOO 120 
SWI00130 
SWIC014 0 
SWI00 150 
SWI00160 
SWI00170 
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c' START OF EACH ITERATION 

/■* -vV -Jr U. Oy* j j- -L- -X> sir sir -W 4r JU -JL- -JL- <JU -dr JL sir JU JU sjlr sb -Jr si, sir X X 

^ i — m ^ t* -i* -t — r^'T-T-'T'-'i- r 

1 LEG= 1 

MTESI=0 
AH=AMASS 
KC OD NT= KCO UN T+ 1 
BORNT=0. 0 

C INITIALIZATION OF STATE (XO) AND COSTATE (QO) 

DO 3 1=1 ,6 
QO (I)=Q01 (I) 

3 X0(I) = X0E(I) 

C INITIALIZATION OF MATRIX OF PARTIALS 

C Z(I,J) PARTIAL OF STATE AND COSTATE WITH RESPECT TO INITIAL COSTATE 
C AND SWITCHING TIMES 

C £(I,J) PARTIAL OF SIGHT END VARIABLES AND SWITCHING CONDITIONS WITH 
C RESPECT TO INITIAL COSTATE AND SWITCHING TIMES 
DO 4 1=1,12 
E(I,13)=0.0 
DO 4 J = 1, 12 
E (I, J) =0.0 

4 Z{I,J)=0.0 
DO 5 1=1,6 

5 Z (1+6,1) =1.0 
LEG5=6 

LEG 6= 7 

IF (ATP (1,4) ) 6,7,6 
C CALLED TO PROPAGATE COAST ARCS 

C PHI (I , J) - PARTIAL OF STATE AT END OF COAST WITH RESPECT TO STATE 
C AT START OF COAST . 

C DPHI ( I, J) - PARTIAL OF COSTATE AT END OF COAST WITH RESPECT TO STATE 
C AT START OF COAST 

7 CALL COAST (PHI , DPHI , OX, LEG, ATP , NO ) 

DO 30 1=1,6 
16=1+6 

DO 30 J=1 , LEG5 
DZ (I, J) =0.0 
DZ ( 16 , J) =0.0 
DO 30 K=1 ,6 

DZ (1, J) =DZ (I,J) +PHI (I ,K) *Z (K , J) 

30 DZ (16 , J) =DZ (16 , J) +PHI (I, K) *Z (K+6,J) +DPBI(I,K) *Z (K,J) 

C UPDATE MATRIX OF PARTIALS Z 

DO 31 1=1, 12 
DO 31 0= 1 , IEG5 

31 Z (I, J) =DZ{I, J) 

0M=DSQRT (QF { 1) **2 + QF (2) ** 2 J -QF (3) **2) 

CBMU=ATP (LE G+ 1 , 4 ) / ( A M * 0 3 ) 

C ADDITION OF NEW COLUMN TC Z CORRESPONDING TO TINS AT END OF COAST 
DO 8 1=1,3 

8 Z (I + 3,LEG6) =-CBMU*QF(I) 

IF (LEG. LE. 1) GO TO 173 

C CALCULATION OF SWITCHING CONDITION AND CORRESPONDING PARTIAL 
E (LEGS , L7) =Ui‘JP-UM 
DO 9 J=1 , LEGS 
E (LEGS , J)=-E (LEG5, J) 

DO 9 K =1,3 


SWI00180 
SWI001 90 
SWI00200 
SWI0021 0 
SWI00220 
SWI00230 
SWI00240 
SWI00250 
SWI 0026 0 
SWIQ0270 
SWI00280 
SNI00290 
SWI 00300 
SWI00310 
SWI00320 
SWI 0033 0 
SWI00340 
SWI 00 35 0 
SWI0036 0 
SWI 00370 
SWI00380 
SWI00390 
SWI00400 
SWI 0041 0 
SWI00420 
SWI00430 
SWI 0044 0 
SWI00450 
SWI0046 0 
SWI0047 0 
SWI 0048 0 
SWI00 490 
SWI00500 
SWI0051 0 
SWI 00520 
SWI00530 
SWI00540 
SWI 0055 0 
SWT00 560 
SWI0057 0 
SWI00580 
SWI00590 
SWI00600 
SWI0061 0 
SWI0062 0 
SPI00630 
SWI00640 
SWI0065 0 
SWI00660 
SWI0067 0 
SWI 00 68 0 
SWI0069 0 
SWI00730 
SWI00710 
SWI 00720 
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9 E (LEG5, J) = E (LEGS, J) + (QF (K) /UJ) *Z (K+6, J) 

173 LEG=L EG+ 1 
LEG5 = LEG+ 5 
LEG6= LEG + 6 
MTEST=1 
GO TO 12 

C CALLED TO PROPAGATE BURN ARCS 
C XO(QO) - STATE (COST ATE) AT START OP BURN 
C XF (QF) - ST ATE (CO ST ATE) AT END OF BURN 
6 CAIL F.KG031 (XO , QO , XF ,QF , Z , EVT , HH»X, LEGMAX,NO) 

BURNT=BURNT+ATP (LS3 + 1 , 1 ) -ATP (LEG, 1 ) 

CALL OUT(XC,QO / XF / QF / LEG'1AX,-1) 

IF (LEG . GE . LSGMAX) GO TO 12 
DH=DSQRT (QF ( 1) *QF(1) + Q? (2) *Q? (2) +QP (3) *QF (3) ) 

UMP=U3 

CBM U=ATP (LEG, 4) / (AH*UM) 

C ADDITION OF NSW COLUMN OP Z CORRESPONDING TO TIME AT END- OF BURN 
DO 14 1=1,3 

14 Z (1+3, LEGS) =CBMUTQF (I) 

C CALCULATION OF PARTIAL OF CONTROL (U) MAGNITUDE WITH RSSP’ECT TO 
C INITIAL COSTATE AND SWITCHING TINES 
DO 15 J= 1 ,LEG6 
DO 15 K=1 , 3 

15 E (L3G6, J) =E (LEG6, J) + (QF (K) /UN) *Z (K+6 , J) 

IF (LEG. LE. 1) GO TO 170 

C CALCULATION OF TR ANSVSR5 ALITY CONDITION AND CORRESPONDING PARTIAL 
C WITH RESPECT TO INITIAL COSTATE AND SWITCHING TIMES 
12 R2=XF (1) ♦XF (1) +XF (2) (2) + <F (3) *XF (3) 

RS=XF ( 1) *QF( 1) +XF (2) *QF (2) + <F (3) *QF (3) 

C3=-UK/ (R2*DSQRT (R2) ) 

C4=-3. 0*C3*SS/R2 

E(L5G5,L7) = XF (U)*Q? (4) + X?(5) *Q? (5)+XF (5) * QF (6 ) -C3*RS 
1 — E (LEGS ,L7) 

DO 16 1=1,3 

DUMMY (I) =— C3*QF (I) -C4 *XF ( I) 

DUMMY (1 + 3) =QF (1+3) 

DUMMY (1+6) =~C3*XF ( I) 

16 DUMMY (1+9) =XF (1+3) 

DO 17 J= 1 , LEG6 

DO 17 K=1 , 1 2 

17 E (LEGS, J) =-DU MMY (K) *Z(K,J) +3 (LEG5,J) 

IF (HTEST. LE. 0) GO T0170 
HTEST=0 

DO 172 J= 1 , LEG6 
172 E (LEGS, J)=-E (LEGS, J) 

E ( L EG5 ,L EG5) =0 . 

LEG=LEG-1 

170 IF (LEG . GE .LEGM AX) GO TO 18 
DO 19 1=1,6 
QO (I)=QF(I) 

XO (I) =XF(I) 

LEG= LEG+1 
LEG5=L 2G+5 
LEG6=LEG J -6 

IF (ATP (LEG, 4) .SQ.O) GO TO 7 


SWI00730 
SWI0Q74 0 
SWI00750 
SWI 00 76 0 
SWI0077 0 
SWI 00780 
SWI00790 
SWI00800 
SWI0081 0 
SWI 00 820 
SWIC0830 
SWT00840 
SWI 00850 
SWI00860 
SW 10087 0 
SWI 0088 0 
SWI00890 
SWI00900 
SWI 009 1 0 
SWI00920 
SWI00930 
SWI00940 
SWI 0095 0 
SW 100 9 60 
SWI 00 970 
SWI 0098 0 
5RIQ0990 
SWI01 000 
SWI01 01 0 
SWIO 1020 
SWI01030 
SWI01040 
SWI01 050 
SWIO 1060 
SPI01070 
SWI01 080 
SWIO 1 090 
SWI01 100 
SWI01 110 
SWI 01 12 0 
SWI01 130 
SWI01 140 
S SI 01 1 5 0 
SWI01 160 
SWI01 170 
SWI01 180 
SWI 01 190 
SWIO 1200 
SPI01210 
SWI01 22 0 
SWIO 1230 
SWIO 1240 
SWI01 250 
SWI01 260 
SWI01 270 
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GO TO 6 

C CALLED TO CALCULATE HI SHI END VARIABLES 
C XF - FINAL STATE , QF - FINAL COSTATE 
18 CALL BVEVAL(XF,QF,Z,C,E,DC) 

DO 20 1=1,6 
16=1+6 

E(I,L7) =DC (I) 

DC (16) =5 (16 , L7) 

20 E (LEG6, 1) =Q01 (I) 

IF(NO.EQ.O) GO TO 79 

C CALLED TO SOLVE SIMULTANEOUS SQJATIONS AND DETERMINE MEN INITIAL 
C COSTATE AND SWITCHING TINES . QO 1 - INITIAL COSTATE 
CALL ADJUST (QO 1 , E , Z , D'J, DUDT , LEGM A X , ATP ) 

CALL OUT (XF, QF, XO , QOI, LEG MAX, 0) 

IF (KCOUNT. GE. KrlAX) GO TO 78 
IF (D'J. GT. .0001 .AMD. DUDT. GT. .001) GO TO 1 
C ALLOWS PROPAGATION 0? NEN TRAJECTORY WITHOUT CALCULATION OF PARTI ALS 
78 NO=Q 

GO TO 1 

79 CALL OUT (XF,Q0E, XO , QF, LEG'l AX, 1 ) 

RETURN 

END 


SRI 01 280 
S RIO 1290 
SWI01 300 
SWI0131 0 
SR 10 1 320 
SWI01330 
SRI01340 
SWI01350 
SRI01360 
S RI0 1 37 0 
SRI01 380 
S RIO 1 390 
SRIO 1400 
SRI01 410 
SWI01420 
SRIO 1430 
SWI01440 
SWI01 450 
SRIO 1460 
SRIO 1470 
SRI 01 480 
SRI01 490 


SUBROUTINE RKG031 (XO,QO, <F , QF ,Z , E VT, H M A X, LEGM AX ,NO) 

IMPLICIT REALMS (A-H,0-Z) 

C THIS ROUTINE GOVERNS TFIS PROPAGATION OF BURN ARCS BY CALLING THE 
C NUMERICAL INTEGRATION ROUTINES. IT DETERMINES THE INTEGRATION STEP 
C SIZE (H) BY AN ESTIMATING ERROR AFTER EACH INTEGRATION STEP. 

C XO AND QO ARE STATE AND COSTA 1 3 AT START OF AFC. XF AND QF ARE STATE 
C AND COSTATE AT END OF ARC. Z IS MATRIX OF PARTIALS. 

COrfMON/CCPIN J/UK , ATP (7,4) , A J ,LEG ,ITSURN 

DIMENSION XO (6) , QO (6 ) , XF (6 ) , QF ( 6) , Z ( 1 2, 1 2) , YN (6 , 13) , YDN (6 , 1 3) , 
1Y3H (6,13) ,YD3fl (6, 1 3) ,Yd(6,13) ,YDM(6,13) ,EY(6) ,EYD(6) 

ITBURN=0 
LEG6=LEG+6 
DO 1 1=1,3 
YN (1,1 )=X0 (I) 

YN ( 1+3 , 1) =Q0 (I) 

YDN (1,1) =X0 (1+3) 

YDN (1+3 , 1) =QO (1+3) 

DO 1 J=2 , IEG6 
J1=J-1 

YN (I, J) =Z (I, J1) 

YN (1+3 , J) =Z (1+ 6, J 1) 

YDN (I, J) = Z (1 + 3, J1) 

YDN (1+3, J) =Z (1 + 9, J1) 

TF = ATP (LEG+1 ,1) 

T0=ATP (LEG, 1) 

H=DSIGN (HMAX ,TF-T0) 

TN=T0 
GO TO 7 
2 E03=H/3. 

XT3URN=ITBURN + 1 


SWI00900 
SWIOO«1 0 
SWI00920 
SWI00930 
SWI00940 
SWI09950 
SWI00960 
SWI00970 
SRI 00 980 
SWI00990 
SUI01000 
SWIC1 01 0 
SWI01020 
SWIP1030 
SWI01 04 0 
SNI01 050 
SRIO 1060 
SWI01070 
SWI01080 
SWI01 090 
SKI01 100 
SWI01 110 
SWI01120 
SWI01 130 
SWI01 140 
SWI01 150 
SWI01 160 
SNI01 170 
SWI01 180 
SWI01 190 
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C A3ALL RKSTEP TO PERFORM NUMERICAL INTEGRATION 
CALL RKSTEP (YN , YDN ,TN, Y3H, YD39 ,H,NO, LEG) 

CALL RKSTEP (YN,YDN, IN, YN , YDM ,H03 ,0 ,LEG ) 

CALL RKSTEP (Y«,YDM,TN+H03 , YM, YDH ,H03 ,0 ,LEG) 

CALL RKSTEP (Yd ,YDd ,TN+HC3*2 YM, YD M, R03 r 0 ,LEG) 

DO 3 1=1,6 

EY (I)=. 1 25E- 1 * (Y8 (I, 1)-Y3ii (1,1) ) 

3 EYD(I) =. 125E- 1- (YD 1 (I,1)-YD3H (1,1)) 

EV2MIN=1 ,E-13*(YDM ( 1 , 1) * YDN ( 1 , 1 } +YDM (2, 1) *YDrt (2,1) + 

1YDM (3 , 1 ) *YDM (3,1) ) 

EVL2=DUAX1 (EVT* (H/ (I? -TO) ) **2,EV2MIN) 

R= (EYD (1}*EYD(1) +EYB (2) * E YD (2) +EYD (3) *EYD(3) ) /EVL2 
C UPDATE STATE , COSTATE AND PARTIALS 
DO S 1=1,6 

YN ( I, 1) =Y 3 (1,1) +EY (I) 

5 YDN (1,1) =YDM (I, 1) + EYD (I) 

DO 6 1=1,6 

DO 6 J=2 , LEG 6 
IN (I,J) =Y3H (I,J) 

6 YDN (I, J) =ZD3H (I, J) 

IF ( DABS (9) .GE. DAB S (TF-TN) ) GO TO 8 
TN=TN +H 

IF (R.LT . 0.04) R = .04 
H=H/R**. 125 

C DETERMINE STEP SIZE (H) FCS NEXT INTEGRATION STEP 

7 IF { DABS(H) .GT. DA3S (TF-TN) ) H=TF— T N 
IF ( DABS (H) .GT. Hd A X) P=DS IGN ( HNAX , H) 

GO TO 2 

8 IF(LEG. IT. LEG 1AX ) GO TO 1C 

C CALCULATION OF PARTIAL OE XF AND QF WITH RESPECT TO FINAL TliJE (TF) 
ON£= 1 .DO 

CALL YDDRHS (YN, YD35, CHS, ATP (LEG + 1,1) ,0) 

L6=6+LEGdAX 
DO 1 1 1 = 1 ,3 
Z (I,L6) =YDN(I, 1) 

Z (1+3 , 16)=YD3H (1,1) 

Z (1+6, 1,6) = YDN (1 + 3, 1) 

11 Z (1 + 9, L6) =YD3H(I+3 ,1) 

10 DO 12 1=1,3 

XF (I) = YN (1,1) 

XF (1 + 3) = YDN (1,1) 

QF (I ) =YN (1+3 , 1 ) 

QF (1+3) =YDN (1+3, 1) 

DO 12 J=2,IEG6 
J1=J-1 

Z (I, J1) =Y N (I,J) 

Z(I + 3, J1)=YDN(I, J) 

Z (1+6, J1) = YN {1 + 3, J) 

12 Z (1+9, J1) =YDN(T+3, J) 

AM=A H— ATP (LEG, 3) * (ATP (LES + 1 , 1) -ATP (LEG , 1 ) ) -ATP (LEG+1 ,2) 

RETURN 

END 


SWI01200 
SNI0121 0 
SWIO 1 220 
SWI01230 
SWI01 240 
SWI01 250 
SWI0 1260 
SWI01270 
SWI01 280 
SWI01290 
SWI01300 
SWI01310 
SWIO 1320 
SWI01330 
SWIO 1340 
SWI01 350 
SWI01 360 
SWI0 1370 
SNI01 380 
SWI01 390 
SWI0 1 40 0 
SWI01410 
SWI01 420 
SWIO 1430 
SWI01 44 0 
SWI01 450 
SWIO 1 46 0 
SWI0 1470 
SWI01 480 
SWI01 490 
SWIO 1500 
SKI01510 
SWI01 520 
SWI01 53 0 
SWI01 540 
SWI01550 
SWI01 560 
SWI01570 
SWIO 1 580 
SWI01590 
SSI01 600 
SWI01 61 0 
SWI01620 
SKI 01 630 
SWI01 640 
SWI01 650 
SWI01660 
SWI01 670 
SWI01 680 
SHI01 690 
SWI01700 
SKI 01 71 0 


A-8 
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SOB ROUTI NE RKST EP { YN , YDN, TN, YN 1 , YDN 1 , H , N, LEG) SWI00010 

IMPLICIT REAL*8 (A-.H, C-Z) SWI00020 

THIS PROGRAM ADVANCES YN AND ID N BY A STEP OF SIZE H TO YN 1 AND YDN1 SWI00030 
USING A FOURTH-ORDER R'JNGF-KUTT A NUMERICAL INTEGRATION SCHEME. IF N SWI00040 
IS POSITIVE, ALL ELEMENTS OF TH'E MATRICES YN AND YDN ARE ADVANCED. SWI00050 
OTHERWISE ONLY THE FIRST COLUMN OF EACH MATRIX IS UPDATED. SWI00060 

DIMENSION YN {6 , 13) , YDN (6,13) ,01 (6,13) ,D2(6,13) , D3 (6,13) ,Y{6,13) SWI00070 

1 , YN1 (6,13) , YDN 1 (6,13) SWI0008Q 

JMAX— 1 SWI00Q9 0 

IF(N.GT.O) JMA X= 6+ LEG SWI00100 

H2= . 5 SWI00110 

CALL YDDRHS (YN,D1 ,H,TN,N) SWI00120 

DO 1 J= 1 , J MAX SWI00130 

DO 1 1=1,6 SWI 0014 0 

1 Y (I, J)=YN (I,J) +H2*(YDN (I,J) +.25*D1 (I,J) ) SWIOO 150 

CALL YDDRHS (Y,D2,H,TN + H2,N) SWI00160 

DO 2 J=1 , JMAX SWI00170 

DO 2 1=1,6 SWI00180 

2 Y (I, J)=YN (I,J) +H* (YDN (I, J) +. 5*D2 (I, J) ) SWI00190 

CALL YDDRHS (Y ,p3, H, T N+H, N) SWI00200 

DO 3 J= 1 , JM AX SWI 0021 0 

DO 3 1=1,6 _ SWI00220 

YN1 (I , J) = YN (I,J) +H* (YDN (I, J) +(D1 (I, J) +2.*D2(I,J))/6.) SWI00230 

3 Y DH 1 (I,J) = YDN (I,J) + (D 1 (I,J) +4.*D2 (I,J) + D3 (I,J) ) /6 . SWI09240 

RETURN SWI00250 

END SWI00260 


SUBROUTINE YDDRHS (Y,YDD,rf,I,N) 

IMPLICIT 3EAL*8 (A-d, C-Z) 

C THIS PROGRAM CALCULATES THE RIGHT HAND SIDES OF THE DIF 
C TIONS DEFINING THE SECOND TIME DERIVATIVES OF R, U AND 
C MATRIX CONSISTING OF RPD, UDD AND WDD . Y CONSISTS OF R, 

C THE STEP SIZE OF PRESENT STEP, T IS PRESENT TIME AND N 
C DETERMINING WHETHER WDD IS TO BE CALCULATED AT THIS EXECUTION 
COMHON/CCPINJ/UA, AT? (7,4) , A 1 ,LEG , ITBUR N 
DIMENSION Y (6,13) ,YDD(6,13) , R (3) , U { 3) ,B (6,6) 

C COMPUTE BASIC QUANTITIES COMMON TO MANY COMPONENTS OF YDD. 
LEG6=6+LEG 
DO 1 1=1,3 
R(I) = Y(I,1) 

1 U(I) =Y(I+3,1) 

AMI =AM— (T- ATP (LEG, 1 ) ) *ATP (LEG, 3) 
H2=1./(R(1)*E(1)+H(2)*B{2) +R (3) *R (3) ) 

02=1. /(U (1) *U (1) +U (2) YU (2) +U (3) *U (3)) 

UM=DSQRT (U2) 

RU=R (1) *U (1) +R{2) * U (2) + 3 (3) * U (3) 

ALPHA=-H*UK*R2*DSQRT (R2) 

BETA=H*uM * ATP (LEG, 4) /A J1 
G.AI!MA=-3 .* ALPHA*R2*R0 
DO 2 1=1,3 

C COMPUTE RDD AND UDD. 

YDD (1,1) =H (I) Y ALPHA + U (I) "BETA 

2 YDD (1 + 3,1) =S (I) YGAMMA+ D (I) *ALPFA 


SWI00270 
SWI 00 280 
FERE NTIAL EQU A-SWI00290 
W. YDD IS A SWT00300 
U AND W, H IS SWI00310 
IS AN INDICATORS FI 00320 


SWI0033 0 
SWIOO 3U0 
SWI00350 
SWI00360 
SWI00370 
SWIOO 380 
SWI00390 
SWI00400 
SWI0041 0 
SWI00420 
SWI00430 
SWI 00440 
SWIOO 450 
SWIOO 460 
SWI00470 
SWI00480 
SKI00490 
SWI 00 500 
SWI 0051 0 
SWIOO 520 
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C 'DECIDE WHETHER HDD IS REQUIRED AT THIS TIME. 

21 IF(N.LE.O) RETURN 

C COMPUTE ADDITIONAL QUANTITIES COMMON TO MANY COMPONENTS OF W DD . 

DE LTA =— 3 . * ALPS A* R2 
EPS IL=— BETA*02 
ZETA=-5. *GAHMA*R2 

C COMPUTE THE MATRIX B NEEDED IN THE MATRIX EQUATION WDD=B*v7 . 

DO 3 J= 1 , 3 
RRJ=DELTA*R (J) 

RUJ=EPSILYU (J) 

URJ=ZETA~2(J) + DELTA A U (J) 

DO 3 1=1,3 

IF (I. EQ . J) GO TO 7 
B {I, J) =3 (!) *RRJ 
B (I, J + 3) =U (I) *3UJ 
B (1 + 3 , J) = H (I) -URJ+U (I) *S8J 
GO TO 3 

7 B (I,J)=R (I) *RRJ + AL?HA 
B (I , J+3 ) =U (I) *R r J J+3 ST A 
B (1+3 , J) =R (I) *UR J+U (I) *RRJ+GA3HA 

3 B (1+3, J+3) =B(I, J) 

C PERFORM THE MATRIX MULTIPLICATION B*W TO GET WDD . 

DO 5 1=1,6 
DO 5 J =2 ,LEG6 
SUM=0 . 

DO 4 K=1 ,6 

4 SUM = SUJ + B(I,K) *Y (K, J) 

5 YDD (I, J) =SUH 

IF (LEG. LE . 1) RETURN 
PRDDd=BETA/AKl 
DO 6 J=8 ,L EG6 

RDDNB=PPDDrt*(ATP(J-7,3) - ATP (J-6 ,3) ) 

DO 6 1=1,3 

C CALCULATE PORTION OF WDD DEPENDING ON DISCRETE CHANGES IN HASS RATE 
6 YDD (I,J)=YDD(I,J) + HD EM 3*0 (I) 

RETURN 

END 


SWI 00530 
SWI00540 
SWI005S0 
SWI 0056 0 
SWI00570 
SWI0058 0 
SWI00590 
SWIOO $00 
SWI00610 
SWI00620 
SWI 0063 0 
SWIQ0640 
SWI00650 
SWI 0066 0 
SWIOO 670 
SWI00680 
SWI 0069 0 
SWIOO 700 
SWI00710 
SWI00720 
SWI00730 
SWIOO 7^0 
SWI 00 750 
SWI 0076 0 
SWI00770 
SWI00780 
SWI 0079 0 
SWI 00800 
SWIOO 810 
SWI00820 
SRI 0083 C 
SWIOO 840 
SRI00850 
SWI0086 0 
SWI00870 
SWI00880 
SWI00890 


SUBROUTINE COAST (PHI, DPHI, UK, LEG , ATP , NO) SWI00010 

IrlPLICIT REALMS (A-H,0-Z) SWI00020 

C THIS ROUTINE PROPAGATES STATE, COSTATE AND THE PARTI ALS OVER COASTING SWI00030 
C ARCS USING A CLOSED FORI UNIVERSAL VARIABLE SOLUTION TO THE TWO-BODY SWI00Q40 
C PROBLEM. PHI = PARTIAL OF ST AT E ( COSTATE) AT END OF ARC WITH RESPECT TOSWI00059 
C INITIAL STATS { COSTATE) . D?HI= PARTIAL OF FINAL COSTATE WITH RESPECT TOSRI0.0060 
C INITIAL STATE. UK IS UNIVERSAL GR AV IAIIONAL CONSTANT. LEG TELLS WHICH SWI00070 
C ARC PROPAGATING. ATP ARRAY CONTAINS SWITCHING TINES AND MASS DATA. SWIC0080 
C IF NO EQUALS ZERO ONLY STATE A AD COSTATS PROPAGATED. SWI00090 

CCIlrfCN/WCO AST/PS Y, ALPHA, FT SWIOO 100 

COMHON/CCOAST/ANN (2,2) ,BNN(3,3) ,XX0(6,6) S4I00110 

1,DAH (2 ,2) ,DBNN (3, 3) ,DXX0 (6,6) - SWI00120 

COHMON/ACOAST/RO (3) ,V0{3) , ERO (3) , DVO ( 3) ,R(3) ,V(3) ,DR(3) ,DV(3) SKI00130 

RO , VO STATE AT START OP COAST , R,V STATE AT END OF COAST SWIOO 1^0 

DR0,DV0 COSTATE AT START OF COAST , DE,DV COSTATE AT END OF COAST SWI00150 


A-1Q 
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DIMENSION PHI (6,6) , DPHI (6,6) , ATP (7,4) ,H0(3) 

C TIME OF COAST 

T=ATP (LEG+1 ,1 ) - ATP (LEG, 1) 

JUKP=0 

RM0=DSQRT (RO (1) *R0 (1) + B0(2) * RO (2) +R0 (3) *R0 (3) ) 

DEM 0= (RO (1) *D80 (1) +R0 (2) *DB9 (2) + R0(3) *DR0 (3) ) /SMO 
SIG0=R0 (1) *V0 (1) + RQ (2) *V0 (2} +30 (3)*V0 (3) 

DSIG0= (VO ( 1 } *DRQ (1) +V0 (2)*DR0(2) +V0 (3) *DPO (3) +R0 (1) *DV0(1) 

1 + 80(2) * DVO ( 2) +50 ( 3 ) *DVO (3) ) 

ALPHA=VO (1) *V0 (1) +V0 (2) -VO (2) +V0 (3) *V 0 (3) -2 ,*UK/RflO 
HO (1) =30 (2) *V0 (3)-R0 (3) *V0 (2) 

HO (2) =R0 (3) *V0 (1)-R0 (1) *V<5 (3) 

HO (3) =S0 (1 ) *V0 (2) -30 (2) *70 ( 1) 

P0= (HO (1) *H0{1) +H0(2) *H0 (2) +H0 (3) *H0 (3) )/0K 
PS Y=T/PO 

ALPSQ=DSQHT (—ALPHA) 

83 ALPSY =PSY*ALPSQ 
SO=DCOS (AIPSY) 

SI =DSIN (ALPSY) /ALPSQ 

C ITERATIVE SOLUTION FOR UNIVERSAL VARIABLE 
S 2= (SO— 1 . 0) / ALPHA 
S3= (S1-PSY) /ALPHA 
FT«RH0*S1 +SIG0*S2+UK*S 3 
RM=3M0*S0+SIG0*S1+UK*S2 
IF (JUMP. £Q. 1 ) GO TO 84 
PSY=PSY+ (T-FT) /SM 

IF (DABS (T-FT) .GE. DABS (T) *1. D-5) GO TO 83 
JU MP= 1 
GO TO 33 

84 FK1=-UK*S2/RI10 

F= 1 . 0+FM1 
FD=-UK*S 1/ (R M*RII0) 

G=FT- UK-S3 
GDh1=-UK*S2/RH 
GD=1 .0+GDM1 
UK R3 = UK / ( R H *R fl *R M ) 

UKR03 = UK/(RMO*RMO*3MO) 

DALPH=2. 0* (VO (1) *D VO ( 1 ) +V0(2)*DV0 (2)+V0 (3) *DV0 (3) +UKR03 * ( 

1R0 (1) *DRO (1) +80 (2) *D SO (2) +R0 (3) *DRO (3) ) ) 

DAP A = DALPH/ALP HA 
DAPA2=DA PA/ALPHA 

DPSY=- (DRKO*S1 + DSIGO*S2+aaO* (PSY*S0-S1) *DAPA* .5+SIGO* (PSY*S1*.5- 
1 52) *DAPA+UK* (PS Y- 1. 5*S 1 +BSY*SO* .5) -DAP A2) /Ml 
DS0=(ALPHA*DPSY+.5*PSY*DAL?H) *S 1 
DS1=S0*DPSY+ (PSY*S9-S1) *DAPA*.5 
DS2=S1*DPSY+ (. 5*PSY*S 1-S2) *DAPA 
DS3=S2*DPSY + (PSY-1 . 5*3 1 + . 5*PSY *S 0) *DAPA2 
S4= (S2-PSY*PSY* .5) /ALPHA 

DSN=S3*DPSY+ ( PS Y *?S Y * » 5-2 . C*S2 + . 5*P SY*S1) *DAPA2 
S5= (S3-PSY*PSY*PSY/6.0)/ALPHA 

DS5=S4*DP5Y+ (PSY*PSY*PSY/6 .0+ (2 . 0*PSY-2 . 5*S 1 + . 5 *PSY *S0) /ALPHA) 
1*DAPA2 

U—S2* ET + UK* (PSY-S4- 3. 0 J -=S5) 

DU=DS2*FT + UK*(DPSY*S4+PSY*DS4-3. 0*DS5) 
DRM=SO*DRMO+DSO*RMO+S1*DSIGO+DS1*SIGO+UK*DS2 


SWIOO 160 
SWI00170 
SWI00180 
SWIOO 190 
SBIOO 200 
S BIO 021 0 
SWIOO 220 
SWIOO 2.30 
SWI00240 
SWIOO 25 0 
SWIOO 260 
SWT0027 0 
SWI00280 
SWIOO 290 
SWI00300 
SWI 0031 0 
SWIOO 320 
SWIOO 330 
SWI0034 0 
SWIOO 350 
SWIOO 360 
SWI0037 0 
SWI 00380 
SWIOO 390 
SWI 00400 
SWI 00410 
SWIOO 420 
SWI00430 
SWI0044 0 
SWI00450 
SWIOO 460 
SWI00470 
SWI00480 
SWIOO 490 
SWI00500 
SWI0051 0 
SWI 0052 0 
SWIOO 530 
SWI 00 54 0 
SWI 00550 
SWIOO 560 
SWI00570 
SWI0058Q 
SWIOO 590 
SWI00600 
SWI0061 0 
SWI0062 0 
SWIOO 630 
SWI 00 640 
SWI 00650 
S KIOO 660 
SWI 00 67 0 
SWI00680 
SWIOO 690 
SWI00700 
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DF= (-UK*DS2-FH1*DBS0) ySMO 

DG=-DK*DS 3 

DG D= (- UK*DS2-GDM1* DRH) /R<1 
EOT=RMO-ea 

DRO 1=RM*DRMO+DP.M*aMO 
D?D={-UK*DS 1- FD*DRO 1) /ROI 
DO 4 1=1,3 

DR (I)=RO(I) *DF + 70(1) *DG+DH0 (I) *F+DV0 (I) *G 
B(I)=RO(I)*F+VO(I)*G 

DV (I) =R0 (I) *DPD+?0 (I) *DGD+DR9 (I) *FD+DV0 (I) *GO 
4 V (I) =HO (I) *FD+ VO (I) *GD 

E,V STATS AT END OF COAST 
DR, DY COSTATE AT END OF COAST 
IF(NO.SQ.O) RETURN 
CALCULATION OF PARTIAL S 

DUKR3-=- 3.0*U KR3*DRM/RM 
DUR03=-3.0*UKR03*D3aO/3dO 
S1RC=S1/RU0 

DS1R0= (DS 1-S1RO*DSMO) /RMO 
S1R=S 1/BM 

DS1 R= (DS1-S1 B* D8M) /Ri 
BO 2=1. 0/ (Rf*0*Rh0) 

B2= 1 .0/ (3M*RH) 

OUK03=- 0*UKR0 3 
DUUK3 = -DU* UK BQ 3-0* UUB0 3 
ANN (1, 1} =-FD*S1 RO-F’l 1*R02 
ANN (1,2) =-FD*S2 
AMN (2,1) =FM1*S1RO+OOK03 
ANN (2, 2) =Fi1 1 *52 
' D£JMM1 = a»N (1,1) . 

DAN (1,1) =-FD*DS 1B0-D?D*S1 R0+OKR03 *DS2+ DUR03*S2 
DU I'M Y= DAN (1 ,1) 

DAN (1,2)=-DFD*S2-?D*DS2 
DAN (2, 1) = Fri 1*DS 1S3+DF*S13C + DUUK3 
DAN (2,2) =FM1*DS2+DF*S2 
CALL A MOLT 
DO 5 1=1,3 
DO 5 J=1,3 
DXXO (I, J) =DBNM (I, J) 

P XXO (I,0)=3NN (X, J) 

IX) 6 1=1,3 

i DXXO (1,1) =DXXO (1,1) +DF 

b XXO (1,1) =XXO(I,I> +F 

ANN (1,1) = 1SN (1,2) 

ANN {2, 1) =ANN (2,2) 

ANN (1,2) =-GDH1*S2 
ANN (2,2) =G*S2-U 
DAN (1 , 1) = DAN (1,2) 

DAN (2, 1) =D AN (2,2) 

DAN (1 ,2} =-GDM1*DS2-DGD*S2 
DAN (2, 2) =-DU+DG*S2+G*DS2 
CALL A MULT 
DO 7 1=1,3 
DO 7 J =1,3 
J3= J+3 


SWI00710 
SWI00720 
SWT00730 
SWI00740 
SWX00 750 
SWI00760 
SWI 00770 
SWI00780 
SWI0079Q 
SWZ00800 
SNI00810 
SWI00820 
SNI00830 
SKI00840 
SWXOO 850 
SKI00860 
SWI00870 
SHI00880 
SWI00890 
SOI 00 900 
SfTIOOSI 0 
SWI00920 
SWI 00 93 0 
SWI00940 
SOI 0095 0 
SWI00960 
SWI00970 
SWI 0^930 
SWI00990 
SWI01000 
SWI01010 
SWIO 1020 
SWIO 1030 
„SWI01 04 0 
SWTO1O50 
SWI01060 
SWI 01 07 0 
SWI 01 080 
SWI01090 
SWI01 100 
SHI01 110 
SWI 01 12 0 
SWI01 130 
SWI01 140 
SWI01 150 
SWI01 160 
SWI01170 
SWI01 180 
SWI01 190 
SWI01200 
SWI01 210 
SWI01220 
SWI 01 23 0 
SWI01249 
SWIO 1250 
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n o 


FILE:’ SRIT<iHD FORTRAN Pi 


CAMBRIDGE MONITOR SYSTEM 


DXXO (I, J3) =DBNN(I, J) 

7 XXO ( 1/ J3) =BNN (I, J) 

DO 8 1=1,3 
13=1+3 

DXXO (I,I3)=DXX0 (1,13) +DG 

8 XX0(I,I3) =XX0 (1,13) +G 
ANN (2,1) =- ANN (1,1) 

ANN (2, 2) =- ANN (1,2) 

ANN (1,1) =-PD*S 1 R-GDM1 *82 
ANN { 1 , 2) = U*UKR 3-GD M 1*S 1 R 
DAN (2, 1) =— DAN (1,1) 

DAN (2,2)=-DAN (1 ,2) 

DAN ( 1 , 1) =-DFD*SlR-?D*DS1R + r JKR3*DS2+DUKR3*S2 
DAN (1,2) =-GDMl *DS1 R-BGD^S 1R +DU-*UKR3+U*DUKR3 
CAIL A MULT 
DO 9 1=1,3 
13=1+3 
DO 9 J=1 , 3 
J3= J+3 

DXXO (I3,J3) =D3NN (I,J) 

9 XX0(I3, J3)=BNN (I,J) 

DO 10 I=li, 6 

DXXO (X , I)= DXXO (1,1) +DGD 

10 XXO (1,1) =XX0 (1,1) +GD 

ANN (1,2) = ANN (1,1) 

ANN (2, 2) =ANN (2, 1) 

ANN (2,1) =-DU MM 1 

ANN (1 , 1) =-FD* (SO/RO 1 +R2 + 3 02) - r JUK03*UKR3 
DAN (1 ,2) =DAN (-1 ,1) 

DAN (2,2) =D AN (2,1) 

DAN (2, 1) =-DU.HY 

DAN (1,1)=-UUK03*DUX33-DaUK3=UXS3-DFD*(S0/R01+R2 + RQ2) -FD* ( (DSO-SO 
1DR01/F01) /R01-2 .0*(R2*DRM/Rd+R02*DRM0/SF0) ) 

CALL A MULT 
DO 11 1=1,3 
X3=I+3 
DO 11 0=1,3 
DXXO (I3,J)=DBNN (1,0) 

11 XXO (13,0) =ENN (1,0) 

DO 12 1=1 , 3 
13=1+3 

DXXO { 13 , L) =DXXO (L3,I) +DFD 

12 XX0(I3,L)=XX0 (13,1) +FD 

XXO PARTIAL OF XF WITH RESPECT TO XO 
DXXO PARTIAL 0? QF 3ITH 33SPECT TO XO 
DO 50 1=1,6 
DO 50 0=1 ,6 
PHI (1,0) =XX0(I,O) 

50 DP HI (I,0)=DXXO (I, J) 

RETURN 

END 

SUBROUTINE AMULT 
IMPLICIT REALMS (A-H,0-Z) 

C SPECIAL PURPOSE MATRIX .HIT PLICATION ROUTINE 
COMMON/CCOAST/ASN (2,2) ,3NN(3,3) ,XXD (6,6) 


SWI01260 
SWIO 1 270 
SWI 01 280 
SWI01290 
SRI01 300 
SWI01310 
SWI01 320 
SWI01330 
SWI01340 
SSI 01 350 
SWIO 1 360 
SWI0 1370 
SWI01380 
S W 1 0 1 39 0 
SWI01400 
SWI01410 
SWI01420 
SWI01 430 
SWI0144Q 
SWI01450 
SWI 01460 
SWI01470 
SWI01480 
SWI01U90 
SWI01500 
SWIO 1510 
SFI01520 
SWI01530 
SWIO 1540 
SWIO 1550 
SWI01560 
* SWI 01 57 0 
SWIO 1 580 
- SWI01590 
SWI 01 600 
SWIO 1610 
SWI01620 
SPI01630 
SWI01 640 
SWIO 1 650 
SWI01660 
SRI 01 67 0 
SWI01 680 
SWI01 690 
SWI01700 
SWJ0171 0 
SRI 01 72 0 
SWIO 1730 
SWI01740 
SWI01750 
SWI 01 76 0 
SWIO 1 770 
SWI01780 
SWI01 790 
SWIO 1800 
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FILEc S»nrCHD FOSTER H Pi 


CAMBRIDGE MONITOR 


SYSTEM 


1 


2 


1, DAN (2, 2) , D3NN (3,3) ,BXX9 (6,6) 

COMKON/ACOAST/RO (3) ,V0 (3) ,DR0 (3) ,DV0 (3) ,E (3) ,V (3) ,DR (3) , DV (3) 
DIMENSION DA (3,2) , A (3, 2] 

DO 1 1=1,3 
DO 1 J=1 , 2 

DA(I,U)=DR (I) *ANN (1 , J) + DV (I) * ANN (2, J) +R (I) *DAN (1, J) +V (I) *DAN{2, J) 
A (I, J) =ANN (1, J) *R (I) +ANN (2,J) *V (I) 

DO 2 1=1,3 
DO 2 J=1 , 3 

DBNN ( J, J) =A (1,1) *DRO { J) + A (1,2) *-DV0 (J) + D A (I , 1 ) +EO ( J) +DA (1,2) -VO (J) 
BNN (I, J) = A{I, 1) 7R0 (J) +A (I, 2) *V0 (J) 

RETURN 

END 


SWI0L81 0 
SWIO 1820 
S WT0 1 ‘830 
SWI01 840 
SWI01 850 
SWI01860 
SWIC1 870 
SWI01 880 
SWI0 1890 
SWI01 900 
SWI01 910 
SWI01920 
SWI01930 


SUBROUTINE ADJUST(Q0 , E , Z, D 0 ,DUDT ,LSG MAX, AT.P) 

IMPLICIT REALMS (A-Ei,0-Z) 

COMrlON/WADJ/A (8) ,CK 

DIMENSION QO (6) ,E(12,13) , Z ( 12, 1 2) , ATP (7, 2) 

C SOLUTION Or SIMULTANEOUS EQUATIONS 
CALL SOLVE (E,LEGMAX) 

L7=LEGMAX+7 

DU=QSQRT (E ( 1 ,L7) { 1 , L7) + E ( 2 ,L7 ) *E (2 , L7) +E (3,L7) *E(3,L7) ) 

DUD=DSQRT(E(4,L7) *E (4,L7) +E (5,L7) *E{5 ,,L7) +E (6,L7)*E (6,L7) ) 
DUDT=DUD* (ATP (LEGMAX+1, 1)-4TP (1, T) ) 

A (1) =. 2/DU 
A (2)'=. 000 3/DUD 
A (3) -1.0 

IF (ATP (1 ,4) .N*E.,0) A (3) = (.5* (ATP (2, 1) -ATP(1, 1) ) /DABS (E (7,L7) ) ) 
CK=QNTN1 (1 .DO, A (1 ), A (2) ,A (3) ) 

C CALCULATION OF BOUND (CK) ON CHANGES ALLOWED IN INITIAL COST ATE 
C AND SWITCHING TIMES 
DO 8 1=2, LEGMAX 
12=1+2 

A (12) = (.5* (ATP (1+1,1)- ATP (1,1) ) /DABS (E (1 + 6 ,L7) - E(I+5,L7) ) ) 

8 CK=DilINl (CK, A (12) ) 

C CALCULATION OF HEW INITIAL COSTATE AND SWITCHING TIMES 
DO 6 1=1,6 

ATP (1+ 1,1) =AT? (1 + 1 , 1) + CK*E (1+6 , L7) 

6 Q0 (I) =Q0 (I) +CK-E (I, L7) 

OM=DSQRT (QO (1) *Q0 (1) +Q0 (2) *Q0 (2) +Q0 (3)*Q0 (3) ) 

DO 7 1=1,6 

7 QO (I) =QC (I) /UM 

RETURN 
END 


SW 101 500 
SWI01510 
SW'IOI 520 
SWIO 1530 
SWI0154 0 
SWI01550 
SWI01 56 0 
SWI0 1 570 
S WI015S0 
SWI 01 59 0 
SWI0 1 60 0 
SWIO 1 6 1 0 
SWI01 620 
SWI01630 
SWI01 640 
SRI01 650 
SWI01 660 
SWIO 1 670 
SWIO 1680 
SWI01690 
SWIO 1 700 
SWI01710 
SWI 01 720 
SKI 01 73 0 
SWI01740 
SWI01750 
SWI 01 76 0 
SWIO 1770 
SWIO 1780 
SWI01790 


SUBROUTINE SOLVE (A, LEG MAX) 
IMPLICIT REAL-8 (A-H,0-Z) 
DIMENSION A (12 , 1 3) 

L6= S+LEGMAX 
L7=7 + LEG MAX 
DO 6 N=1,L6 


SWI01800 
SWI01 81 0 
SWI 01820 
SWIO 1830 
SWI01840 
SWI01 85 0 
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FIL E« SWITCHD FORTRAN Pi 


CAMBRIDGE 


MONITOR SYSTEM 


IBIG=N 
DO 1 I=H,L6 

IF ( DABS ( A (I, N ) ).GT. DABS (A (IB IG, H) ) ) IBIG=I 

1 CONTINUE 

IF (IBIG.EQ.N) GO TO 3 
DO 2 J=N,L7 
Q= A (N, J) 

A(ff, J) —A ( IBIG , J) 

2 A(I3IG,J)=Q 
3 DO 5 I— 1 ,L6 

IF (I. SO- N) GO TO 5 
Q=A(I,N)/A(N,N) 

M=N+1 

DO 4 K=3,L7 

4 A (I r K) =A (I,K) -Q*A (N , K ) 

5 CONTINUE 

6 CONTINUE 

DO 7 1*1,16 

7 A ( I, L7) =A ( I, L7) /A (I f I) 

RETURN 

END 


SWI 01 860 
SWI01870 
SWI01880 
SWI01 890 
SWI01 900 
SWI01910 
SWI01920 
SWI01930 
SWI 01 94 0 
SWI01950 
SWI01960 
SWI 01 970 
SWI01980 
S NIC 1990 
SWI 02 000 
SWI 02 01 0 
SWI02020 
SWI02030 
SWI02 040 
SWI 02 05 0 
S WI02060 


C 

r* 

C 

c 

c 

c 

r 

n 


SUBHOUIINE BVEVAL(XF,QF,Z,C,E, DC) 8VE0Q010 

IUPLICIT REAL*8 (A-H,0~Z) BVE00020 

BVEQ0030 

THIS VERSION DEALS WITH A FIVE-CONSTRAINT RIGHT-END BOUNDARY- VA LUE BVE0004Q 

PROBLEM WHERE THE FIVE CONSTRAINED FUNCTIONS ARE THE THREE COMPON- BVE00050 
ENTS OF THE ORBITAL ANGULAR VELOCITY VECTOR (R CROSS V) AND THE FIRSTS VE0006P 
TWO COMPONENTS OF THE VECTOR /(HOSE DIRECTION IS THE DIRECTION OF BV500070 

PSPICENTER AND WHOSE *1 A GNITUD3 IS THE ORBITAL ECCENTRICITY. THUS, INBVE00080 
EFFECT, ALL OF THE SIX CLASSICAL ORBITAL ELEMENTS ARE CONSTRAINED SVE00090 
EXCEPT THE MEAN ANOMALY, WHICH IS FREE. ' BVE00100 

BVE0011 0 


COMHOK/CCPINJ/OK, ATP (7,4) , A I ,LEG , ITBUR N 
DIMENSION XF <6 ) ,C ( 12) ,DC (1 2) , G (7, 6} , Z ( 12 , 1 2) ,E (1 2, 1 3) ,R (3) , V (3) 
1,QF(5) , DUMMY ( 12) 

LEG 5= 6+ LEG 
DO 1 1=1,3 
R(I> = XF ( I) 

1 V(I)=XF(I+3) 

R2=R (1) **2 + R ( 2) *#2 + R (3) **2 
G(1,2)=V(3) 

G (1 ,3) =- V (2 ) 

G (2,3) =V ( 1 ) 

G(1 ,5) =-R (3) 

G(1,6) =B{2) 

G (2,6) =-R (1) 

RM=DSQST (R2) 

R3=RM*S2 

C1=-1 . Q/RA + (V ( 1) ** 2 + V (2) *-*2+V (3) **2) /UK 
C2=- (R (1) *V (1) +R (2) *V (2) + R (3) *V (3) )/UK 
DO 4 1=1,3 
DO 3 J= 1,3 


BVEOO 120 
BV50013Q 
BVE00140 
BVEOO 150 
BVEOO 160 
BVE00170 
BVE0018Q 
BVEOO 190 
BVEO0200 
BVEQ021 0 
BVEOO 220 
BVE00230 
BVE00240 
BVE00250 
BVEOO 260 
BVE00270 
BVE00280 
BVEOO 290 
BVE00300 
BVEC031 0 
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FILE: • SKITGHD FORTRAN Pi 


CAMBRIDGE MONITOR SYSTEM 


I? (I. LE. J) GO TO 2 
G(I,J) =-G(J,I) 

G (I, J+3) = -G (J,T+3) 

2 G {1 + 3, J) =R (I) *R (J) /R 3— V (I) *V (J) /UK 

3 G (1+3, 0 + 3) = (S (I) *V (J) -2 . -V (I) *R (J) } /UK 
G (I,I)=0. 

G (1,1+3) =0. 

G (1+3 , 1} =G (1+3,1) +C1 

4 G(I+3,I+3)=G(I+3,I+3) +C2 
DO 5 1=1 ,3 

DC (I)=0. 

DO 5 J= 1 , 3 

5 DC(I)=DC(I) +G(I,J)*3(J) 

DO 8 1=1, 2 
SUH=0. 

DO 7 J=1 ,3 

7 SUM=SUM + G (I,J)*DC{J) 

8 DC (1+3) =C (1+3) +R (I) /BM+ SOM/UK 
DO 9 1=1,3 

9 DC(I) =C(I) -DC (I) 

DO 10 1=1,5 
DO 10 J=1 , LEG6 
E(I, J) =0. 

DO 10 K=1 , 6 

10 E (I,J)=E (I, J) +G (I,K) *Z(K,J) 

SS=XF( 1) (1) +XF (2) *QF (2) + XF (3) *QF(3) 

C3=- UK/R3 
C 4=-3. 0*C3-RS/32 

DC (6)= XF (4) *QF (4) +XF(5) # QF (5) +XF(6) *QF (6) -C3*-RS 

DO 16 1=1,3 

DUMMY (I)=-C3*QT(I) -C4*XF(I) 

DUMMY (1 + 3) =QF (1+3) 

DUMM Y ( 1 + 6) =- C 3 X F (I) 

16 DUMMY (1 + 9) =XF (1 + 3) 

DO 17 J=1,LSG6 

DO 17 K= 1 , 1 2 

17 E (6, J) =-DUMMY (<) *Z (K,J) + E {6, J) 

EETURN 

END 


BVEOO 320 

BVE00330 

BVE 00340 

BVEOO 350 

BVE00360 

BVE00370 

BVEOO 380 

BVEOO 39 0 

BVE00400 

BVE004 10 

BVE00420 

BVE00430 

BVE00440 

BVEOO 450 

BVEOO 460 

8VE00470 

BVE00480 

BVEOO 490 

BVE00500 

BVE0051 0 

BVEOO 520 

BVE00530 

BVE00540 

BVE0055 0 

BVEOO 560 

BVE00570 

BVE00580 

BVEOO 590 

BVE00600 

BVE0061 0 

BVEOO 620 

BVE00630 

BVE00640 

BVE0065 0 

BVEOO 660 

BVE00670 

BVE00680 

BVE00690 

BVEOO 700 


SUBROUTINE BVEVAL (XF ,Q F ,Z , C, E, DC) 

IMPLICIT REAL-8 (A-H,0-Z) 

THIS VERSION OF BVEVAL IS FOR 6 CONSTRAINT RENDEZVOUS MISSIONS 
THE FINAL ORBIT IS SPECIFIED BY A RADIUS , VELOCITY AND TIME 
CONTAINED AS INPUT CONSTANTS IN THE C VECTOR 

COMNON/CCPINJ/DK ,ATP (7,4) , AN, LEG ,1 TBURN 
COMOON/O IN/PL ANS (54) ,NOUT, LOGIC 

DIMENSION XF{6) ,QF (6) r Z (12 ,12) ,C (1 2 ) , E ( 12 , 1 3) , DC (1 2 ) , RO (3 ) , VO (3) 
1,H0 (3) ,V (3) ,R (3) , E 1 (3) 

PROPAGATE RADIUS AND VELOCITY OF FINAL 03 BIT TO FINAL TIME 
T=ATP (LSG+1, 1) -C (7) 


BVEOO 01 0 
BVE00020 
RVE00030 
BVE0004 0 
BVE0005C 
BVE00060 
BVE00070 
BVE00080 
BVE00090 
BVEOO 100 
BVEOO 1 1 0 
BVE001 2 0 
BVEOO 130 
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PILE.; SvJITCHD FORTRAN 


PI 


CAMBRIDGE MONITOR SYSTEM 


DO 40 1=1,3 
RO (I)=C (I) 

40 70 (I) =C (1+3) 

JUilP=0 

RM0=DSQHT ( RO ( 1) *R 3 ( 1) + RO ( 2) * RO (2) +R0 (3) *R0 (3) ) 

SIG0=R0 (1) *V0 (1) +50 (2) *70(2) +50 (3) *70 (3) 

ALPHA=V0 (1) *V0 (1) +V0 (2) *70 (2) +70 (3) *70 (3) -2.*UK/PM0 
HO (1) =50 (2) *7 0 (3) -5 0 (3) *70 (2) 

HO (2) =R0 (3) *70 (l)-SO (1) *70 (3) 

HO (3) =50 (1) *70 (2) -5 0 (2) *70 (1) 

P0= (HO (1) *H0 (1) +IIO (2) *H0 (2) +H0 (3) *H0 (3) ) /UK 
PS Y=T/PO 

ALPSQ = DSQRT (-ALPHA) 

83 ALPSY =PSY*ALPSQ 
SO=DCOS (ALPSY) 

Sl=DSIN (ALPSY) /ALPSQ 

C ITERATIVE SOLUTION FOR UNIVERSAL VARIABLE 
S2=(S0-1.0) / ALPHA 
S3= (S1-PSY) /ALPHA 
FT=RMO*S 1 +SIGO *S 2+ U K*5 3 
RE=RM0*S0+SIG9*S1+UK*S2 
IF (JOHP.EQ . 1) GO TO 34 
PS Y=PSY+ (T-FT)/RM 

IF ( DABS (T- FT) .GE.DABS (T)*1 .D-5) GO TO 83 
JUMP=1 
GO TO 83 
84 CONTINUE 

Fnl =— UK*S2/EM0 
F=1 . 0+FH1 
FD=— U K*S 1/ (Rtf*RM0) 

G=FT- UK *5 3 
GDM1--UK*S2/RM 
GD= 1.0 +GD Ml 
DO 4 1=1,3 

R (I) =R0 (I) *F+ VO (I) *G 
V (I) =S0 (I) *FD+ VO (I) *GD 
DC ( I) =R ( I) — XF (I) 

4 DC (1 + 3) =V (I)-XF (1+3) 

R2T=R ( 1) **2+R ( 2) **2 + S(3) **2 
S3 T=H2T* D5 QRT (R2T) 

L FG5=L EG+5 
DO 5 1=1,6 
DO 5 J= 1 , LEGS 

5 E (I, J) = Z (I, J) 

DO 6 1=1,3 

E (I, LEG5 + 1 ) = Z (I , LEG5 + 1 ) -V (I) 

6 E (1+3, LEG 5+1) =Z (1+3, LEG 5+1) +5 (I) *UK/R3T 
WRITE (NOUT, 100) R,V 

100 FORM AT (20X, 42HTARGST VEHICLE FINAL POSITION AND VELOCITY, /25X , 
1 2HR , 3E 1 4. 6, 5X, 2HV ,3214.6) 

IF {LOGIC. GT. 0 ) RETURN 
LOGIC= 1 

El (1) «- (BO (1) /RH0+ (HO (2) *V0 (3)-H0 (3) *V0 (2) ) /UK) 

El (2) =- (RO (2) /5.U0+ (HO (3) *V0 (1) -HO (1) *V0 (3) ) /UK) 

El (3) =- (SO (3) /RH0+ (HO (1) *V0 (2)-H0 (2)*V0 (1 ) ) /UK) 


BVE00140 

B VEOO 150 

BVE00160 

BVE00170 

BVEOO 180 

BVE00190 

BVE00200 

BVE00?10 

BVE00220 

BVE0023 0 

BVE00240 

BVE00250 

BVE00260 

BVE00270 

BVEOO 280 

BVE00290 

BVE00300 

BVE0031 0 

BVEOO 320 

BVE00330 

BVE 0034 0 

BVEOO 350 

BVE00360 

BVE0037 0 

BVEOO 380 

BVE00390 

BVE00400 

BVEOOA1 0 

BVEOO 420 

BVE00430 

BVE0044 0 

BVEOO 450 

BVE0046 0 

BVE00470 

BVE00480 

BVEOO 490 

BVS00500 

BVE0051 0 

BVEOO 520 

BVE00530 

BVE00540 

BVE 0055 0 

BVEOO 560 

BVE00570 

BVE00580 

BVE 00590 

BVEOO 600 

BVE0061 0 

BVE00620 

BVEOO 630 

BVE00640 

BVE00650 

BVE 0066 0 

BVE0067 0 

BVE00680 
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TILE; SvTITCHD FORTH AS ?1 


CAMBRIDGE MONITOR SYSTEM, 


E«AG=DSQRT {El (1) **2+El (2) **2+E1 (3) **2) 
fill AG=DSQRT{P0*TfK) 

AAXIS=-UK/ ALPHA 
ENERGY=. 5*ALPHA 

PER 10 D= (6. 2831853) *D5QRT (DABS (A AXIS**3/UK) ) 

RhIN=AAXIS* (1.0-EMG) 

RM AX=A AXIS* (1 . O+EMAG) 

WRITE (NOUT ,101) AAXI3 , RMI N, EM AX, ENER GY ,P EP IOD , HM A G, HO , EMAG , E 1 
101 FORMAT ( 22X r 27 HI AR GS T ORBIT SPECIFICATIONS ,/25X, 15 HSEMIMAJOR AXIS 
IE 14. 6,1 X5HRHIN= f E14.6 , 1 X53BnA«X= ,E14. 6, 1 X7 HEN33GY= , El u . 6 ,/25X, 

27 HPER IOD = , E 1 A . 6 , 1 X5H rfdAG=,El4.6 ,1 X8HH VECTOR, 3E1 4 . 6 / /25X5HE?1AG=, 
3E14 .6,1 X8HE VECTOR, 3E14.5) 

BETORN 

END 


BVE00690 
BVE00700 
B VE 0071 0 
BVE00720 
BVE00730 
BVE00740 
BVE00750 
BVE00760 
, BVE0077 0 
BVEC0780 
BVE0079 0 
BVE00800 
8VE0081 0 
B VS 0082 0 
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APPENDIX II-l 

INPUT FILE and OUTPUT LISTING 
for 5-CONSTRAINT MISSION 



FILE: FILE 


FT 02 FOCI ° l 


CAMBRIDGE MONITOR SYSTFM 


1 


TEST CASF NUMBER 2 ro* DOUBLE PRFriSION DECK 

B-fONSTR AINT , 44 OcGRFE 3 LANE CHANGE CASE RUN f)N 7/7/69 
APRIL 28,1970 

398601.5 12644651. 1Q0. 

4551 .3085947 19 .8398425. 05 763245. 59906 10 -5. 41 70895- . 01 1 838° 

.4601788 -.8868545 .0413273 -.0004394 -.0010504 -.0004081 

65248.406 62230.797 -93156.688 


410 

934. 0.0 

2047.6868 0.0 
23 02.86 77 0.0 
21032. 055 0.0 
21160.352 0.0 


0.0 0.0 

22384.406 92936.438 

0.0 0.0 


22384.406 92006.438 

0.0 0.0 



CASE NUMBER 1 


TEST CASE NUMBER 2 FOR DOUBLE PRECISION DECK 

5-CONSTRA INT t 44 DEGREE PLANE CHANGE CASE RUN ON 7/7/69 
APRIL 28 t 1970 

GRAVITATIONAL CONST ANT= 398601.50 INITIAL MASS= 0.12644651D 08 MAX INTEGRATION STEP* 100.000 
MAX NUMBER OF IT£RATIONS= 10 NUMBER OF SEPERATE ARCS= 4 


ATP ARRAY 


TIME 


DELTA MASS 

MASS RATE 


THRUST 


0.93400000D 

03 

0.0 

0.0 


0.0 


0 .204768680 

04 

0*0 

0.223644060 

05 

0.92986438D 

05 

0.230286770 

04 

0*0 

0.0 


0*0 


0.210320550 

05 

0.0 

0 *223844060 

05 

0*929864380 

05 

0.211508520 

05 

0.0 

0.0 


0.0 



INITIAL STATE XO 
ESTIMATED 00 


4551.1085900 

0.4601788 


4719.83 98400 
-0.8868545 


25. 05763Z4 
0.0415273 


5.5990610 

-0.0004394 


-5.4170895 -0.0113389 

-0.0010504 -0.0004081 


DESIRED FINAL C 0.652484060 05 0.62230797D 05 -0.931566880 05 0.0 


0.0 


0.0 


0.0 


ITERATION NUMBER 1 

COAST ARC 

LEG= 1 STATE AT END 0.5S6318D 04 -0.3267670 04 -0.356783D 01 -0.391783D 01 -0.6759350 01 -0. 3192020-01 

COSTATE AT END -0.722406D 00 -0.129057D 01 -0.3238770 00 -0.1143170-02 0.648873D-03 -0.1473300-03 

PSY- 0.170086D 00 ALPHA= -0.60^8930 02 CALCULATED COAST TIME* 0.11t369D 04 

SEMIMAJOR A XI S= 0.654634D 04 RMIN= 0.6531210 04 RMAX= 0.6561460 04 ENERGY= -0.3044463 02 
PER 100= 0.527118D 04 HMAG= 0.510819D 05 H VECTOR 0.798617D 02 0.194I82D 03 -O.510BI5O 05 

EMAG= 0. 2310290-02 E VECTOR 0. 8254970-04 -0.2308790-02 -0.8647610-05 RMAG= 0.6538290 04 
BURN ARC " ' 

LEG= 2 STATE AT END 0.4263510 04 -0.5045860 04 -0.745324D 02 -0.7175290 01 -0.7217930 01 -0.5838540 00 

COSTATE AT END -0.969136D 00 -0 .107409D 1)1 -0.3460460 00 -0.7569490-03 0. 1007890- 02 -0.2598350-04 

MASS AT END OF LEG= 0.693258D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST ARC 

LEG= 3 STATE AT END -0.3176180 05 0.2596260 05 0.30269LD 02 0.111210D 01 0.1199750 01 0.9415060-01 

COSTATE AT END -0.6258110-01 -0.3221320 00 0.1580300 01 0. 2782580-04 -0.6460410-05 0. 1008000-04 

PSY= 0.761753D 00 ALPHA* -0.167 481D 02 CALCULATED COAST TIME= 0.187292D 05-- 

SEMIMAJOR AXIS= 0.2379980 05 RMI N= 0.654657D 04 RMAX= 0.4105200 05 ENERGY* -0.8374073 01 
PERIOD* 3.3654000 05 HMAG= 0.670907D 05 H VECTOR 0.2408080 04 0.302406D 04 -0.6697920 05 
EMAG= 0.724931D 00 E VECTOR 0.571933D 00 -0.445440D 00 0.4512110-03 RMAG= 0.4102280 05 

BURN ARC \ 

LEG= 4 STATE AT END -0.316326D 05 0.2608210 05 0.1489523 03 0.105819D 01 0.7803870 00 0.206199D31 

COSTATE AT END -0.5925830-01 -0.3228990 00 0.158143D 01 0.2811410-04 -0.644141D-05 0.8994333-05 

MASS AT END OF LEG= 0. 4273380 07 NUMBFR OF INTEGRATION STEPS IN BURN ARC* 2 

RESULTANT ORBIT SPECIFICATIONS 

SFMIM4J0R AXI S= 0.296053D 05 RMIN= 0.176885D 05 RMAX= 0.415220D 38 FNERGY= -0.6731933 01 
PE RI00= 0 . 506949 D 05 HMAG= 0.994420D 05 H VECTOR 0.536648D 05 0.653337D 05 -0.522855H 05 
EMAG= 0.4025210 00 E VECTOR 0.3309460 00 -0.219748D 00 0.6487880-01 RM4G= 0.4099900 05 

TOTAL BURN T I ME* 0.373978D 03 ARC TIMES 0.111369D 04 0.25518ID 03 0.1872920 *5 0.1187970 33 

DC 0.1158360 05 -0.3152930 04 -0.4087120 05 -0.3309460 03 0.219748D OO 3.676^670-05 

-0 . 29 L28 30-03 -0.126509D 00 -0.1915530-04 0.0 0.0 O.o 

DETERMINANT OF E= 0.469171D 06 DIAGONAL OF E 0.343838D 06 3.5835570 06 0.^544090 04 -0.7473240 03 

-0.4484330 03 0.3764130 03 -0. I 150340-01 0.2915420-04 -0 .440405 D-03 0.2773570-07 

DU* 0.7131250-01 DUDT= 0.208593D 01 CK= 0.1000000 01 EVT= 0 . 1 300000-37 
r K = M T N *V" O 4 <s,0 01 n .Z n ">'"0 '',rif) / '3 *>l r i , v-ir> <3? *.270670 01 n. 520883 0* 



CHANGE REQUESTED IN INITIAL COSTATE, SW ITCHING TINES AND FINAL TIME 

-0.1 33239D-0 1 -0. 1268250-01 -0.677398D-01 -0.5533B 10-04 -0 . 4391 040-04 -0. 7483 1 1 0-0 A 
0.7850230 02 0.709078D 02 -0.3388870 04 -0.337746D 04 

END OF ITERATION NUMBER 1 

NEW 00 0 .44073 ID 00 -0.897259D 00 -0. 261461D-01 -0. 4939840-03 -0. 1091540-02 -0.4817080-03 

NEW SWITCH TIMES J. 2126190 04 0.237378D 04 0. 176432D 05 0. 177734D 05 

ITERATION NUMBER 2 

COAST ARC 

LEG= l STATE AT END 0.533120D 04 -0.378316D 04 -0.604647D 01 -0.453386D 01 -0.636427D 01 -0 . 31281 5D-01 

COSTATE AT END -0 .9923850 00 -0.129064D 01 -0. 405562D 00 -0. 130251 D-02 0.8372060-03 -0.4187680-04 

PSY= 0.1 82 094 D 00 ALPHA= -0.608893D 02 CALCULATED COAST TIME” 0.1192190 04 

SEMIMAJOR AX I S= 0.6546340 04 RMI N= 0.653I21D 04 RMAX= 0.656146D 04 ENERGY” -0.3044463 02 
PERIOD” 0.5271180 04 HMAG= 0.510819D 05 H VECTOR 0.798617D 02 0.194182D 03 -0.510815D 05 

EMAG= 0. 2310290-02 E VECTOR 0.8254970-04 -0.230879D-02 -0.8647610-05 RM AG= 0.6537130 04 
BURN ARC 

LEG” 2 STATE AT END 0.3819450 04 -0.536248D 04 -0.776274D 02 -0.7806200 01 -0.6400160 01 -'>.5957890 30 

COSTATE AT END -0. 12691 7D 01 -0.102912D 01 -0.398199D 00 -0 . 8982 16D-03 0.123767D-02 0.1083133-03 

MASS AT END OF LEG” 0.710258D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC” 3 

COAST ARC 

LEG” 3 STATE AT END -0.252695D 05 0.2535680 05 0.737217D 02 0.131941D 01 0.129997D 01 0.1ini84D 00 

COSTATE AT END 0.229815D 00 0. 297048D 00 0.1707950 01 9 . 355 534D-0 4 -0.393605D-04 -3 . 7595530-05 

PSY= 3. 708845 D 00 ALPHA” -0.188264D 02 CALCULATED COAST TIME” 0. 1526940 35 

SEMIMAJOR AXIS” 0.2117250 05 RMI N= 0.654643D 04 RMAX= 0.3579860 05 ENERGY” -0.9413199 01 
PERIOD” 0.3065970 05 HMAG” 0.664230D 05 H VECTOR 0.269808D 04 0.288156D 04 -0.663055D 05 

EMAG” 0.6908050 00 E VECTOR 0.488843D 00 -0.488100D 00 -0. 132043D-02 RMAG= 0.357983D 05 
BURN ARC 

LEG” 4 STATE AT END -0.250786D 05 0. 2554620 05 0.215366D 03 0.163953D 01 0.164O57D 01 0.2253330 01 

COSTATE AT END 0.234419D 00 0.291910D DO 0.170683D 01 0.3516700-04 -0.395685D- 04 -0. 952655D-05 

MASS AT END OF LEG” 0.418797D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC” 2 

RESULTANT ORBIT SPECIFICATIONS 

SEMIMAJOR AXIS” 0. 337465D 05 RMIN” 0.316603D 05 RMAX” 0.358327D 05 ENERGY” -0.5905823 01 
PERIOD” 0.6169560 05 HMAG” 0.115758D 06 H VECTOR 0.572109D 05 0.568635D 05 -0.8302690 35 

EMAG” 0.618197 D-0 1 E VECTOR 0.373544D-01 -0.4866900-01 -0.7592820-02 RMAG” 0.3579930 05 
TOTAL 8 URN TIME” 0.377793D 03 ARC TIMES 0.119Z19D 04 0.247586D 03 0.152694D 05 0.130207D 03 

DC 0.8037550 04 0.536732D 04 -0.1012980 05 -0.373544D-0 1 0. 4866900-01 -0. 11818 ID-04 

0.1098620-04 -0 . 669537D-01 -0. 2280 43D- 04 0.0 0.0 0.0 

DETERMINANT OF E= 0.904181D 06 DIAGONAL OF E 0. 165&63D 06 -0.2682310 06 -0.77010OD 03 -0.5708990 03 
0. 348 1330 03 -0.2007700 04 0.143213D-01 0.2719790-04 0.3834920-02 0. 4432910-07 
DU” 0.4647370-01 DUDT= 0. 186750D 01 CK= 0.10O000D 01 EVT” 0.100000D-07 
C K=MI N OF 0.43035D 01 0.27051D 01 0.100Q00 01 0. 141910 02 0.22269D 01 0.12645D 02 

CHANGE REQUESTED IN INITIAL COSTATE, SWITCHING TIMES AND FINAL TIME 
-0.2437840-01 -0.1306300-01 0.3734770-01 0.275731D-04 -0 . 1654470-04 0.1061360-03 

0 . I66599D 02 0.253832D 02 0. 3453808 04 0.3448660 04 

END OF ITERATION NUMBER 2 

NEW QO 0.4159030 00 -0.9093400 00 0.IH895D-01 -0.465903 0-03 -8 .1186890-82 -0.375167 0-0 3 

NEW SWITCH TIMES 0.214285D 04 0.239916D 04 0.2109700 05 0.2122200 05 



ITERATION NUMBER 3 


COAST ARC 

LEG= 1 STATE AT END 0.525461D 04 -0.388844D 04 -n.656638D 01 -0.465°69D 01 -0.6273080 01 -0.3113160-31 

COSTATE AT END -0.104533D 01 -0.130990D 01 -0.312801D 00 -0. 130499D-02 0.879293D-03 -0 . 627670D-04 

PSY= 0. 184643D 00 ALPHA= -0.608893D 02 CALCULATED COAST TIME= 0.120885D 04 

SEMIMAJOR AXI S= 0.654634D 04 RMIN= 0.653121D 04 RMAX= 0.656146D 04 ENFRGY= -0.304446D 0? 

PERIOO= 0 .8271 18 D 04 HMAG= 0.510819D 05 H VECTOR 0.798617D 02 0.194182D 03 -0.5108150 05 

EMAG = 0. 2 3L029D— 02 E VECTOR 0 .825497D-04 -0 .2308790-0 2 -0 .8647610-05 RMAG= 0.6536390 04 

BURN ARC 

LEG= 2 STATE AT END 0.364010D 04 -0.5495030 04 -0.6719150 02 -0.8081350 01 -0.6273090 01 -0.484941D 00 

COSTATE AT END -0.1327610 01 -0.102640D 01 -0.314121D 00 -0. 860 159D-03 0.128918D-02 0.519496D-04 

MASS AT END OF L EG= 0.6907310 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST ARC 

LEG= 3 STATE AT END -0.305264D 05 0.3042630 05 0.2605320 02 0.102638D 01 0.1179740 0] 0.7*73820-01 

COSTATE AT END 0.4243350 00 0.379562D 00 0.158398D 01 0.393361D-04 -0.3555950-04 0.5514150-05 

PSY= 0.751374D 00 AIPHA= -0.160457D 02 CALCULATED COAST TIME= 0.1869780 05 

SEMIMAJOR AXI S= 0.248416D 05 RMIN= 0.654756D 04 RMAX= 0.4313570 05 ENERGY= -0.8022860 01 
PERIOO= 0.389655D 05 HMAG= 0.6731900 05 H VECTOR 0.2243270 04 0.230823D 04 -0.672420D 05 

EMAG= 0.736428D 00 E VECTOR 0.5088200 00 -0.532379D 00 -0.1300290-02 RMAG= 0.431QO1D 05 
BURN ARC 

LEG= 4 STATE AT END -0.3036570 05 0.3060030 05 0.151467D 03 0.1592810 01 0.1644300 01 0.2106040 01 

COSTATE AT END 0.429241D 00 0.375098D 00 0.1584610 01 0.391126D-04 -0. 3583650-04 0.4528030-05 

MASS AT END OF L EG= 0.4107950 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 2 

RESULTANT ORBIT SPECIF ICAT IONS 

SEMIMAJOR AXIS* 0;452129D 05 RMI N= 0.429755D 05 RMAX= 0.4745020 05 ENEPGY= -0.4408050 01 
PERIOD* 0 .9567620 05 HMAG= 0.134081D 06 H VECTOR 0.6419650 05 0.641927D 05 -0.986707D 05 

EMAG= 0. 4948430-01 E VECTOR -0.4182490-01 0.2365650-01 -0.1182150-01 RMAG= 0.4311010 05 

TOTAL BURN TIME* 0.38L368D 03 ARC TIMES 0.1208850 04 0.256310D 03 0.186978D 05 0.1250530 03 

DC 0.105187D 04 -0.196191D 04 0.551401D 04 0.4182490-01 -0.2365650-01 0.6361640-05 

-0.7953360-05 0.2406340-01 0. 143150D-04 0.0 0.0 0 

DETERMINANT OF E= 0.2188250 07 DIAGONAL OF E 0.9894280 05 0.104159D 07 -0.4048480 03 -0.8459240 03 

-0.3021590 04 -0.1730550 04 -0.2217750-01 0.1247790-04 -0.119960D-02 0.2922160-07 

0U= 0.219740D-02 DUOT= 0.4744450 00 CK= 0.100000001 EVT= 0.100000D-07 

CK=MIN OF 0.91017D 02 0.128280 02 0.10000D 01 0. 13542D 03 0.57158D 03 0.421190 03 

CHANGE REQUESTED IN INITIAL COSTATE, SWITCHI NG T1 MF S AND FINAL TIME 
-0.61 2783D-03 -0 . 305946D-03 -0.2087930-02 -0 . 193894D-05 -0.3072180-05 -0.2310160-04 
0.3653560 01 0.270719D OX -0.1364910 02 -0.1379750 02 

, END OF ITERATION NUMBER 3 

NEW QO 0. 4152900 00 -0.9096440 00 O'. 9 10 1570-02 -0. 467846D-03 -0. 1109960-02 -0. 3982670-03 

NEW SWITCH TIMES 0.214650D 04 0.240187D 04 0.2108330 05 0.2120820 05 



teratjon m umber 4 


COAST ARC 

LEG= 1 STATE AT END 0*5237540 04 -0*3911320 04 -0*6680060 01 -0.4687040 01 -0*6252750 01 - 
COSTATE AT END -O.105682D 01 -0.131006D 01 -0 .332533D 00 -0 . 1308Q6D-0 2 0*8881900-03 - 
ni . M PSY= °«165201D 00 ALP HA= -0*608893D 02 CALCULATED COAST TI ME= 0.1212500 74 

SEMIMAJOR AKIS= 0.6546340 04 RMIM= 0.6531210 04 RMAX- 0.656146D 04 ENFROY= -0. 3044460 02 
£ ERKI0 = 0-527119D 04 WMG= 0. 5108190 05 H VECTOR 0.798617 0 02 0.1941820 03 -0.51^5 0 05 

EMAG= 0. 2310290-02 E VECTOR 0.825497D-04 -0.2308790-02 -0.8647610-05 RMAG= 0.6536840 04 


31 0971 D- 01 
61 53970-04 


BURN ARC 

LEG= 2 STATE AT END 
COSTATE AT END 


0. 362387D 04 -0.550437D 04 -0.69T0980 02 -0.8092320 01 
-0.133891D 01 -0.1025370 01 -0.3326820 00 -0.8630440-03 


•0.6232320 01 -0. 
0 . 1 29790D-02 0. 


COAST ARC 
LEG= 


KM AC f A T . „ “ . _ ~-A.MAA.~~,,,, ~ a K, . . c KA C. KM UU ~KJ A O O 3 t U“U 1 U.I 

MASS AT END OF LEG- 0.692849D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 3 


3 STATE AT ENO 
COSTATE AT END 


-0.2909330 05 
0.362727D 00 


0.3047240 05 
0.3432090 00 


0.6882270 02 0.1141140 01 0.H1211D 01 0. 

ncv- n i cimnn nr, ~ ; ” — 0.1644490 0 1 0.337 839D-04 -0. 349777D-04 0. 

ccu.u^Tn™ A 2rl 56l n 2 n 00 ALPHA= -0.163768D 02 CALCULATED COAST TIME= 0.186815D 05 

0.243394D 05 RMIN= 0.6547350 04 RMAX= 0.4213150 05 ENERGY= -r. 8188430 pj 
PERIOD= 0.3778990 05 HMAG= 0.672126D 05 H VECTOR 0.2356870 04 0.2401820 04 -0.671283D 05 

BURN ARC °* 730998D °° E VECTQR 0- 5027790 00 -0.5306320 00 -0 . 1333730-02 RMAG= 0.42^060 Ss 

LEG= ^ nclirc A J !^ D -0*2892360 05 0.306344D 05 0.1960810 03 0.1616380 01 0.151T500 01 0. 

OO 0*3388360 00 0.1644450 01 0*3354050-04 ~0 .352D42D-0* -0. 

RESULTAN? OR81T l"lcT ? IC Tint' ^ ^ INTEGRA7I0N STEP S ™ «*" AR0 = * 

n!l?J^ JaR ^ A A IS= °* 42l672 ° 05 RMIN= 0.4210570 05 P,MAX= 0.4222870 05 ENERGY= -0.4726449 01 
c«A^- D= n ?/'ltiIn 2 5o°£ . HMAG= 0* 1296450 06 H VECTOR 0.650717D 05 0.6203560 05 -0.934083D 05 

TnTM E nnou -r?u=^ 8 rt 6 oIno-»f V ECT0R -0 * 1204020-02 0.244341D-04 -0.8225350-03 RMAG= 0.4213170 05 

TOTAL BURN 7INE= 0.380273D 03 ARC TIMES 0.1212500 04 0.2553630 03 0.186815D 05 0.1249100 03 


5071120 OO 
507607 ( 9-04 


7985640-01 

2351470-06 


2133850 01 
8596B7D-06 


° C n o °* 1952000 03 0.2516540 03 0.120402D-02 -0.2443410-04 -0.9670360-06 

-9.2287010-06 0.3041250-03 -0.3383360-06 0.0 0.0 0.0 

S! °I DIAG3NAL 0F e 0.7064140 05 0.102992D 07 -0.1037600 04 0.5794690 03 

0.230050D 04 -0.391689D 03 -0.2095080-01 -0 . 103926D-04 -0.8341490-02 0.2790160-07 

. 562 91 0-02 0U0T= 0. 1908 270-01 CK= O.IOOOOOD 01 EV7= 0. 1 OOOOOD-07 

'CK MIN 3F- 0 ,10'1890— 03' 0.318730-03“ 0. 10000D 01 0.76645D 04 0.199530 03 0.654270 03 


]END OF ITERATION NUMBER 


CHANGE REQUESTED IN INITIAL C OST ATE ,S W ITCH! NG T IMPS AND FINAL TIMF 
-0.3020010-03-0.1185070-03 0.1935920-02 0.854805D-06 0.3912910-06 

-0 .9511650 00 -0.903245D 00 0.459106D 02 0.4600610 02 


0.459566D-07 


NEW GO 0.4149870 00 -0.9097600 
NEW SWITCH TIMES 0.214555D 04 


00 0.1103750-01 -0. 4669900-03 -0.1109560-02 - 0 . 3987200-03 

0.240096D 04 0.2U293D 05 0.212543D 05 



ITERATION NUMBER 5 


COAST ARC 

L£G= 1 STATE AT END 0.5241990 04 -0.3905370 04 -0. 6659480 01 -0.4679930 01 -0.6258060 01 -9.3119610-91 

COSTATE AT END -0.1054760 01 -0.13I074D 01 -0.3321880 00 -0.1308860-02 0.8861030-0? -0. 6A?740D-O4 
PSY= 9. 185056 D 00 ALPHA= -0.608893D 0? CALCULATED COAST TIME= 0.121155D 34 

SEMIMAJOR AX1 $= 0.654634D 04 RMIN= 0.653121D 04 RMAX= 0.656146D 04 ENERGY® -0.3044463 02 

PERIOD® 0.5271180 04 HMAG= 0.510819D 05 H VECTOR 0.798617D 02 0.1941820 03 -0.5108150 05 

EMAG= 0.231029D-02 E VECTOR 0.8254970-04 -0,2308790-02 -0 . 864761 D-05 RMAG= 0.653635D 04 


LEG= 2 STATE AT END 0.3629780 04 -0.550056D 04 -0.6971940 02 -0.808573D 01 -0.62M79D 01 -0.5075350 OO 

COSTATE AT END -0.1336950 01 -0.1026500 01 -0.3330380 00 -0.8634310-03 0.1296080-02 0.5704800-04 

MASS AT END OF LEG® 0.692 7420 07 NUMBER OF INTEGRATION STEPS IN BURN ARC® 3 

COAST ARC 

LEG® 3 STATE AT END -O t 29I240D 05 0.304933D 05 0.705079D 02 0.1142960 01 0.110834D 01 9. 7984410-91 

COSTATE AT END 0.358020D 00 0.339462D 00 0.164571D 01 0.3344250-04 -9.3470960-04 0.4356750-06 

PSY® 0. 757596D 00 ALPHA® -0. L63647D 02 CALCULATED COAST TIME® 9.187283D 05 

SEMIMAJOR AXIS® 0. 243574D 05 RM IN® 0.654736D 04 RMAX® 0.4216740 05 ENERGY® -0.8182350 
PERIOD® 0.3783180 05 HMAG® 0.6721650 05 H VECTOR 0.2356560 04 0.2405970 04 -0.671321D 95 

EMAG= 0.731196D 00 E VECTOR 0.503535D 00 -0.530187D 00 -0. 1 325790-02 RMAG® 0.4216T0D 05 


LEG® 4 STATE AT END -0 . 289 543D 05 0.306547D 05 0.198116D 03 0.161297D 01 0.1599740 01 0.2138520 01 
COSTATE AT END 0 .362185D 00 0.335109D OO 0.1645693 01 0. 3320 20D-04 -0. 349324D-04 -9. S 5899 9 D-9 6 

MASS AT ENO OF LEG® 0.412925D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC® 2 

RESULTANT ORBIT SPECIFICATIONS 

SEMIMAJOR AXIS® 0.421738D 05 RMIN® 0.421675D 05 RMAX® 0.421802D 05 ENERGY® -0.4725709 01 
PERIOD® 0. 861936D 05 HMAG® 0.129656D 06 H VECTOR 0.652565D 05 0.6223870 05 -0.931584D 05 

EMAG® 0 . 1503 5 8 D-03 E VECTOR -0.1102900-03 0.1017730-03 - 0.9262770-05 RMAG® 0.4216750 05 
TOTAL BURN TIMS® 0.380417D 03 ARC TIMES 0.121I55D 04 9.255411D 03 0.1872830 05 0.125005D 93 


DC -0.8115180 01 — 0 . 794244D 01 0.167323D 01 0. 1 102900-03 -0 . 1 01 7730-03 0.271221D-08 
9. 2571740-08 0.8532200-04 0.140462D-09 0.0 0.0 0.0 

DETERMINANT OF E= 0.1745620 07 DIAGONAL OF E 0.703419D 05 0.104080D 07 -0.102216D 04 0.567671D 93 
-0.2988640 04 -0.3187600 03 -0.2065410-01 -0 . 122545D-04 -0 . 88318 ID-02 0.276103D-07 
DU= 0.7954620-05 OUDT® 0. 275393D-P3 CK= 0.1000000 01 EVT= 0.100000D-07 
C<=MIN OF 0 .25 143 D 05 9.221360 05 0.10000D 01 0. 191050 06 0.26600D 05 9.475140 04 


CHANGE REQUESTED IN INITIAL COSTATE, SWITCHING TIMES AND FINAL TIME 

0.7298420-05 0.3296250-05 0.670297D-06 -0 . 73 1 178D-08 0.509422D-08 -0. 1 921 080-07 

—0.2 16644D-02 -0.1498000-02 0.3505320 00 0.337377D 00 

END OF ITERATION NUMBER 5 

NEW 00 0.4149940 09 -0.909757D 00 0.110381D-01 -0 . 466997 D-03 -0.1109560-02 -0. 3°82310-03 

NEW SWITCH TIMES 0.2145550 04 0.240096D 04 0.211296D 05 0.2125460 05 



I TERATION NUMBER 6 


COAST ARC 

LEG= l STATE AT END 0.5242000 04 -0.390535D 04 -0.6650410 01 -0.4679910 01 -0.625837D 01 -0.3] lOfelD-01 

COSTATE AT END -0.105475D 01 -0.131074D 01 -0.332196D 00 -0.1308850-02 0.8860960-03 -0. 64277 ID-04 

PSY= 0.1950560 00 ALPHA= -0.608893D 02 CALCULATED COAST T IME= 0.1211550 34 

SEMIMAJOR AX IS= 0.654634D 04 RM IN= 0.653121D 04 RMA X= 0.656146D 04 FNERGY= -0.3044463 ■ 02 
PE R 100= 0.527118D 04 HMAG= 0.510819D 05 H VECTOR 0.798617D 02 0.1941820 03 -0.51O815D 05 

EMAG= 0.231029D-02 E VECTOR 0.825497D-04 -0. 230879D-02 -0. 8647610-05 RMAG= 0.653635D 04 
BURN ARC 

LEG= 2 STATE AT END 0.362978D 04 -0.550056D 04 -0.697216D 02 -0.808573D 01 -0.624172D 01 -0.5075530 00 

COSTATE AT END -0.133694D 01 -0. L02649D 01 -0.333047D 00 -0. 863 42 8)3—03 0.129607D-02 0.570483D-04 

MASS AT END OF L EG= 0.692741D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST ARC 

LEG= 3 STATE AT END -0.291242D 05 0.304936D 05 0.705177D 02 0.114297D 01 0. 11^8320 01 n.7984610-01 

COSTATE AT END 0.35R071D 00 0. 3394820 00 0.164577D 01 0.334445D-04 - 0. 34711 ID- 04 0.048240-06 

P$Y= 0. 757602D 00 ALPHA= -0.163646D 02 CALCULATED COAST TIME= 0.L87286D 35 

SEMIMAJOR AX IS= 0. 2435760 05 RMIN= 0.654736D 04 RMA X= 0.421678D 05 ENERGY= -0.8182283 01 
PERIOD^ 0.378323D 05 HMAG= 0.6721650 05 H VECTOR 0.235664D 04 0.240606D 04 -0.671321D 05 
EMAG= 0.73119BD 00 E VECTOR 0.5035380 00 -0.5301880 00 -0 . 132578D-02 RMAG= 0.421674D 0^ 

BURN ARC 

LEG= 4 STATE AT END -0.289545D 05 0.306550D 05 0.198097D 03 0.161296D 01 0.150966D 01 0.213823D 31 

COSTATE AT END 0.362236D 00 0.335129D 00 0.164575D 01 0.332039D-04 -0. 349339D- 04 -0. 6587450-06 

MASS AT END OF LEG= 0.412953D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 2 

RESULTANT 0R81T SPECIFICATIONS 

SEMIMAJOR AXI S= 0.421680D 05 RMIN= 0.421679D 05 RMAX= 0.421680D 05 FNERGY= -0.4726353 01 
PERIOD= 0.361756 D 05 HMAG= 0.129647D 06 H VECTOR 0.652485D 05 0.6223090 05 -3.9315670 05 
EMAG= 0.122783D-05 E VECTOR -0 .8 396160-06 0.895822D-06 0.1034910-07 RMAG= 0.^216790 05 

TOTAL BURN TIME= 0.380404D 03 ARC TIMES 0.121155D 04 0.255412D 03 0.187286D 05 0.124992D 03 



TEST CASE NUMBER 2 FOR DOUBLE PRECISION DECK 

5-CONSTRAINT , 44 DEGREE PLANE CHANGE CASE RUN ON 7/7/69 
APRIL 28 ,1970 


SUMMARY TABLES 


ITERATION TOTAL BURN TIME LENGTH OF BURN AND CD AST ARCS 

NUMBER 


1 

0. 37397790 

03 

0. 1 1 13687D 

04 

0.2551809D 

03 

0.1872919D 

05 

0.11 87970D 

03 

2 

0 .37779350 

03 

0 • i 192 1890 

04 

0.2475864D 

03 

0.1526941D 

05 

0 • 1302070D 

03 

3 

0.3813682 D 

03 

0 • 1208849 D 

04 

0 .25630980 

03 

0. 186978 3D 

05 

0 .12505840 

03 

4 

0.38027330 

03 

0. 1 2 125 03D 

04 

0. 2 553634D 

03 

0.1868147D 

05 

0 • 1249099D 

Q3 

5 

0.38041670 

03 

0.12115510 

04 

0.25541130 

03 

0. 1872829D 

05 

0. 1250054D 

03 

6 

0 . 3804042D 

03 

0. 12115490 

04 

0.2554120D 

03 

0.1872864D 

05 

0*1 249922D 

03 


ITERATION 

NUMBER 


ERROR IN BOUNDARY CONDITIONS-DCU-8 > 


1 

2 

3 

4 

5 


-0.115836D 
0.803 755D 
0.105187D 
- 0.1767030 
-0.8115L8D 


05 -0,3152930- 04 -0.408712D 
04 0. 536732 D 04 -0.101298D 

04 -:0. 19619,10 04 0.55]L401D 

03 9.-19520^0-03 - 0.2516540 

01 -0.7942440 01 0.167323D 


05 -0.330946D 00 0.219748D 00 0.676267D-05 -0 .291 283D-03 -0.1265090 00 
05 — 0.373544D-0 1 0 .486690D-01 f -0. 1 1 818 10-04 0.10 98620-64,-0.6695370-01 
04 0. 4182490-01 -0.236565D-Q1 Q^,636164D-D5 -0 .795336D-05 0. 240634D-01 
03 0.1204020-02 -0.2443410-04 -0 ; 567036D-06 -0. 228701D-06 0.3041256-03 
01 0. I L0290D— 03 -0. 101773D-03 0.2712210-08 0.257174D-08 0. 8532200-04 


ITERATION 

NUMBER 


1 

2 

3 

4 

5 




d'ai9r-i2) , 

V* 

00 

DUDT 

CK 


0. 19L5531D— 04 

0.0 

a*o 

- 0.0 

0.71312 54D- 01 . 

0.20859340 01 

0* 100Q3000 

01 

-0. 2280426D— 0^ 

0. 0 

d.o ' 

o.d 

0 .46473670-01 

0.186749 8D 01 

0. 10000003 

01 

0.14315010-04 

0.0 

0.0 

0.0 

0. 21974000-02 

0. 47444530 00 

0.10000003 

01 

-0. 33833 56D-06 

0.0 

0.0 

d.o 

0 * 19629 13D-02 

0. 19082710-01 

0. 10000003 

01 

0.14046250-09 

0.0 

0.0 

0.0 

0 .79546 16 D- 05 

0.2753928D-03 

0.10000003 

01 


l 

ITERATION NEW QO GENERATED BY PRESENT ITERATION 

‘NUMBER 

l 


0 

1 

2 

3 

4 
s 


0 .46017 88 OD 00 
0 *440730710 00 
0 *415903370 00 
0*415289590 00 
0.414986790 00 
0.4 1499400D 00 


-0.886854500 00 
-0.397258500 00 
— 0 *909339970 00 
-0.909643730 00 
-0.909760480 00 
— 0 • 9097571 9D 00 


0 *415273000—01 
-0.261461310-01 
0.111895240-01 
0 .9 10 15700D-02 
0.110374670-01 
0.110381370-01 


-0. 43940 0000-03 
-0.493983660-03 
-0.465907710-03 
-0. 467845520-03 
-0 « 46698 98 2 P-0 3 
-0 • 466 9971 3D-03 


-0. I p5 040000—02 
-0.109 15385D-02 
-0.110688850-02 
-0.110995000-02 
-0 . 1 1095 6460-0? 
-0. I \ 0955 95D— 0? 


-0.408100000-03 
-0. 481707830-03 
— 0 . 37516651 D-03 
-0.398267110-0* 
-0 .398220390-03 
- 0 . 308230600-^3 



APPENDIX I I- 2 

INPUT FILE and OUTPUT LISTING 
for RENDEZVOUS MISSION 



FILE: FILE 


FT02F0 0 1 PI 


CAMBRIDGE MONITOR SYSTEM 


PAGE 001 


1 

TEST CASE BOMBER 2 FOR DOUBLE PRECISION DECK 
RENDEZVOUS , 44 DEGREE PLANE CHANGE (7/8/69) 

APRIL 28,1970 
398601.5 12644651. 100. 

4551.3088 4719.843 25.057641 5.5990612 -5. 41 70902- . 01 18389 

.4601788 -.8868545 .0415273 -.0004394 -.0010504 -.0004081 

-28954.51430655.047 198.07783 1.6129526 1.5096618 2.1382224 21254.526 
410 

934. 0.0 0.0 0.0 

2047.6868 0.0 22384.406 92986.438 

2302.8677 0.0 0.0 0.0 

21032.055 0.0 22384.406 92986.438 

21150.852 0.0 0.0 0.0 



CASE NUMBER 1 


TEST CASE NUMBER 2 FOR DOUBLE PHECISION DECK 
RENDEZVOUS , 44 DEGREE PLANE CHANGE (7/8/69) 

APRIL 28,1970 

GRAVITATIONAL CONSTANT= 398601. 5C INITIAL MASS 3 0.12644651D 08 MAX INTEGRATION STEP= 100.000 
MAX NUMBER OF ITERATIONS 3 10 NUMBER OF SEPERATE ARCS 3 4 


ATP ARRAY 


TIME 


DELTA MA S3 

MASS RATE 


THRUST 


0.93400000D 

03 

0.0 

0.0 


0.0 


0.20476868D 

04 

0.0 

0*223 84406D 

05 

0.9 2986438D 

05 

0.23028677D 

04 

0.0 

0.0 


0.0 


0.21032055D 

05 

0.0 

0 . 2238440 6D 

05 

0.929B64380 

05 

0 . 21 150852D 

05 

0.0 

0.0 


0.0 



INITIAL STATE XO 
ESTIMATED QO 


4551.3088000 

0.4601788 


471 9; 8430000 
-0.8868545 


25.0576410 5.5990612 

0.0415273 -0.0004394 


-5.4170902 
-0. 0010504 


-0.0118389 
-0.000408 1 


DESIRED FINAL C -0. 28954514D 05 0.30655047D 05 0.19807783D 03 0.16129526D 01 0.15096618D 01 0.21382224D 01 0.21254526D 05 


ITERATION NUMBER 1 

j» " O ' 4 Jt5 ' ( % 

COAST ARC *«. 

LEG 3 1 STATE AT END 0.566318D 04 -0.326767D 04 -0.356782D 01 -0.391782D 01 -0.675935D 01 -0. 318202D-01 

COSTATE AT END -0.722407D 00 -0.129057D 01 -0.323877D 00 -0 . 1 14 317D-02 0.648871D-03 -0. 147330D-03 

^PSY- 0.170086D 00 ALPHA 3 -0.608892D 02 CALCULATED COAST TIME 3 0.111369D 04 

,\ SEMI MAJOR AXIS 3 0. 654634D 04 RMIN 3 0.653122D 04 RHAX 3 0.656147D 04 ENERGY 3 -0.304446D 02 
' PERIOD^ 0.527119D 04 HBAG 3 0.510820D 05 H VECTOR 0.798618D 02' 0.194182D 03 -0.510815D 05 
EH AG 3 0.231018D-02 B VECTOR 0.831756D-04 -0 .230866D-02 -0 .86461 3D-05 RMAG 3 0.653829D 04 
BURN ARC 

LEG 3 2 STATE AT END 0.426352D 04 -0.504586D 04 -0.745324D 02 -0.717528D 01 -0. 721794D 01 -0.589854D 00 

COSTATE AT END -0.969137D 00 -0.107410D 01 -0.346046D 00 -0 . 756952D-03 0.100789D-02 -0.2598380-04 

MASS AT END OF LEG 3 0.693258D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC 3 3 

COAST ARC 

LEG= 3 STATE AT END -0.317621D 05 0.259625D 05 0.302587D 02 0.111208D 01 0.119976D 01 0.941503D-01 

COSTATE AT END -0 . 6 26 227D- 0 1 -0.322131D 00 0.1580300 01 0.278217D-04 -0. 645883D-0S 0 100B99D-04 

PSY= 0.761752D 00 ALPHA 3 -0.167481D 02 CALCULATED COAST TIME 3 0.1B7292D 06 

SEMIMAJOR AXIS 3 0.237998D 05 RMIN 3 0.654658D 04 RMAX= 0.410531D 05 ENERGY 3 -0.8374040 01 
PERIOD 3 0.365402D 05 BMAG= 0.670907D 05 H VECTOR 0.240808D 04 0.302406D 04 -0.669793D 05 

EMAG= 0. 724932D 00 E VECTOR 0.571934D 00 -0.445440D 00 0.451239D-03 R MAG 3 0.410230D 05 

BURN ARC 

LEG= 4 STATE AT END -0.316329D 05 0.260821D 05 0. 148941D 03 0.10S811D 01 0.780402D no 0.206199D M 

COSTATE AT END -0 . 59300 5D- 0 1 -0.322897D 00 0.158143D 01 0.281100D-04 -0.643479D-05 0 899516D-05 

MASS AT END OF LEG 3 0.427338D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 2 

TARGET VEHICLE FINAL POSITION AND VELOCITY 

R -0.291209D 05 0.304977D 05 -0.236038D 02 V 0.159695D 01 0.152651D 01 0.213827P 01 

TARGET ORBIT SPECIFICATIONS 

SEMIMAJOR AXIS 3 0. 4216780 05 RMIN 3 0.421677D 05 RMA X= 0.421679D 05 ENERGY 3 -0.472637D 01 
PERIOD 3 0.86 1752D 05 H MAG 3 0.129646D 06 H VECTOR 0 .652483D 05 O,622307D 05 -0.9115670 4^ 

EMAG= 0.266008D-05 E VECTOR 0. 170 160D-05 -0 . 203762D-05 -0 . 169346D-06 
RESULTANT ORBIT SPECIFICATIONS 

SEMIMAJOR AXIS 3 0.296051D 05 RMIN 3 0.176883D 05 RM AX= 0.415220D 05 ENERGY 3 -0.673196D 01 
PERIOD 3 0.506945D 05 HMAG= 0.994416D 05 H VECTOR 0.536647D 05 0.653842D 05 -0.522841 D ns 
EMAG= 0.402527D 00 B VECTOR 0.334948D 00 -0.219758D 00 0.648665P-01 RMAG 3 0 . 4099923 05 

TOTAL BURN TIME 3 0.373978D 03 ARC TIMFS 0.1 11369D 04 0.255181D 03 0 18779?D 95 0 1 137Q7D 4-) 



DC 0.251 195D 04 0 . 44 1 559D 04 -0. 17 2545D 0 3 0.538835D00 0.746112D00 0.762818D-01 

-0. 29 1278D-03 -0.126512D 00 0.191551D-04 0.0 0.0 0.0 

DETERMINANT OF E® -0.492687D 10 DIAGONAL OF E -0.252892D 05 -0.160285D 04 0.803543D 02 0.59O142D 08 

-0.176823D 05 -0.245039D 04 -0.222588D 00 -0.333908D-04 -0 . 2 1 249 3D-0 1 0.3745700-05 

00= 0.101254D 00 DUDT= 0.165794D 01 CK= 0.1000000 01 FVT= 0.100000D-07 

CK=MIN OF 0.19752D 01 0. 36S82D 01 0.10000D 01 0.69923D 03 0.25975D 03 0.119230 o 2 

CHANGE REQUESTED IN INITIAL COSTATE, SWITCHING TIMES AND FINAL TIHE 
— 0. 5655 14D-0 1 -0. 32961 ID-01 -0.7725 14D-0 1 -0. 166398D-04 -0 .713886D-04 -0. 367702D-04 
0.873701D 02 0.875525D 02 0.123604D 03 0. 118622D f>3 

END OF ITERATION NUMBER 1 

NEW QO 0.401574D 00 -0.915137D 00 -0.355424D-01 -0. 453720D-03 -0. 1 1 1608D-02 -0.442607D-O3 
NEW SWITCH TIMES 0.213506D 04 0.239042D 04 0.211557D 05 0.212695D 05 

ITERATION NUMBER 2 

COAST ARC 

LEG® 1 STATE AT END 0.529070D 04 -0.383939D 04 -0.632351D 01 -0.460105D 01 -0.631605D 01 -0. 31 2033P-01 

COSTATE AT END -0. 106080D 01 -0.133603D 01 -0.374985D 00 -0. 134022D-02 0.870025D-03 -0.204363D-04 

PSY® 0. 183450D 00 ALPHA= -0.608892D 02 CALCULATED COAST TIME® 0. 120106D 04 

SBMIMAJOB AXIS® 0.654634D 04 RMIN= 0.653122D 04 RMAX® 0.656147D 04 ENERGY® -0.3044460 02 
PERIOD® 0.527119D 04 HMAG= 0.510820D 05 H VECTOR 0.798618D 02 0.194182D 03 -0.510815D 05 

EMAG= 0.23101 8D-02 E VECTOR 0.831756D-04 -0. 230866D-02 -0.864613D-05 R MAG® 0.653700P 04 
BURN ARC 

LEG® 2 STATE AT END 0.369857D 04 -0.545243D 04 -0.746042D 02 -0.800957D 01 -0.632648D 01 -0.849005D 00 

COSTATE AT END -0.135158D 01 -0.105384D 01 -0.362839D 00 -0 . 897869D-03 0. 129670D-02 0.114184D-D3 

MASS AT END OF LEG® 0.692849D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC® 3 

COAST ARC 

LEG® 3 STATE AT END -0.287717D 05 0.302228D 05 0.102612D 03 0.118885D 01 0.108231D 01 0.871027D-01 

COSTATE AT END 0.261560D 00 0.275796D 00 0.179783D 01 0.291824D-04 -0. 325724D-0U -0. 107309D-04 

PSY® 0.761312D 00 ALPHA® -0.165123D 02 CALCULATED COAST TIME® 0.187652D 05 

SBMIMAJOB AXIS® 0.241397D 05 RMIN® 0.654735D 04 RMAX® 0.417320D 05 ENERGY® -0.825615D 01 
PERIOD® 0.3732560 05 HMAG® 0.671694D 05 H VECTOR 0.252143D 04 0.262808D 04 -0.670705D 05 
EHAG® 0 . 728772D 00 E VECTOR 0.506813D 00 -0.523686D 00 -0 . 146705 D-02 RM AG® 0.417282D 05 
BURN ARC 

LEG® 4 STATE AT END -0.286210D 05 0.303599D 05 0.210473D 03 0.147978D 01 0.13&725D01 0 195011D 01 

COSTATE AT END 0.264870D 00 0.272081D 00 0.179655D 01 0.289874D-04 -0. 32711 ID-04 -0. 1 1 8530D-04 

MASS AT END OP LEG® 0.438081D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC® 2 

TARGET VEHICLE PINAL POSITION AND VELOCITY 

R -0. 289304D 05 0.306776D 05 0.230041D 03 VO 161525D 01 0. 150722D 01 0.213821D 01 

RESULTANT ORBIT SPECIFICATIONS 

SSMIMAJOR AXIS® 0.352788D 05 RMIN® 0.288256D 05 RMAX® 0.417319D 05 ENERGY® -0.5649310 01 
PERIOD® 0.659448D 05 HMAG® 0.116583D 06 H VECTOR 0.589215D 05 0.561255D 05 -0.834855D 05 
EMAG= 0.182918D 00 E VECTOR 0.129191D 00 -0.129428D 00 0.416719D-02 RM AG® 0. 417245D 05 

TOTAL BORN TIHE® 0.369178D 03 ARC TIMES 0.120106D 04 0.255363D 03 0.187652D 05 0.113815D 03 

DC -0 . 309362 D 03 0.317683D 03 0. 195678D 02 0. 135471D 00 0. 159980D 00 0.1 88098D 00 
0. 176786D-04 -0. 857154D-01 -0. 22446 ID-04 0 0 0.0 0.0 

DETERMINANT OF E= 0.345157D 10 DIAGONAL OF E -0.207878D 05 -0.14S524D 04 0.618760D 02 -0.850652D 05 
-0.824959D 04 -0.225935D 04 -0.300208D 00 0.144148D-03 0.528743D-02 0.508284D-05 
DU® 0.555362D-01 D UDT= 0.102966D 01 CK= 0.100000D 01 EVT® 0.100000D-07 
CK=M IN OF 0.36013D 01 0.59249D 01 0.10000D 0 1 0. 12304D 05 0.241 80D 03 0.50<>47U 01 

CHANGE REQUESTED IN INITIAL COSTATE ,S WTTCHI NG TIMES AND FINAL TIME 

0. 177390D-0 1 0 . 57524 0D-02 0.523 1 17D-0 1 -0. 1 76409D-04 0.774860D-05 0.468245D-04 

0.107927D 02 0.108031D 02 -0.280006D 02 -0.168308D 02 

END OF ITERATION NUMBER 2 

NEW QO 0.418668D 00 -0.907985D 00 0.167435D-01 -0. 470636D-03 -0.1106630-^2 -0 . 3 9S 17U0-0 3 

NEW SWITCH TIMES 0.214585D 04 0.240122D 04 0.211277D 05 0.212526D 05 



ITERATION NUHBER 3 


COAST ARC 

LEG 3 1 STATE AT END 0.524060D 04 -0.390723D Oa -0.665974D 01 -0.468215D 01 -0 6256400 01 -0 311933D-M 

COSTATE AT END -0.104918D 01 -0.130487D 01 -0.328919D 00 -0 . 1 30 332D-0 2 0.8R511UD-03 -0 . 704827D-0 4 

PSY 3 0.18S101D OO ALPHA 3 -0,60889 2D 02 CALCULATED COAST TINE 3 0,1211850 04 

SBHIHAJ08 AXIS 3 0.654634D 04 RNIN= 0.653122D 04 RMAX= 0.656147D 04 ENERGY 3 -0.3044460 02 
PERIOD 3 0.527119D 04 HSAG= 0.510820D 05 3 VECTOR 0.798618D 02 0.194182D 03 -0.510815D 05 

ENAG= 0.231018D-02 E VECTOR 0.831756D-04 -0. 230866D-02 -0. 864613D-06 RNAG= 0.653685D 04 
BURN ARC 

LEG 3 2 STATE AT END 0.362823D 04 -0.550177D 04 -0.695469D 02 -0.808656D 01 -0.6239911) 01 -0.50651011 00 

COSTATE AT END -0.133009D 01 -0.102130D 01 -0.331484D 00 -0 . 859594D-03 0.129238D-02 0.4993430-04 

BASS AT END OF LEG= 0.692826D 07 NUHBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST ARC 

LBG= 3 STATE AT END -0.290958D 05 0.304904D 05 0.704041D 02 0.114524D 01 0.110708D 01 0. 7971940-01 

COSTATE AT END 0.392962D 00 0.353330D 00 0.163520D 01 0.354282D-04 -0. 3610300-04 0 1430250-05 

PSY = 0.757655D 00 ALPHA 3 -0.163720D 02 CALCULATED COAST TIME= 0. 187264D 35 

SEHIM AJOE AXIS 3 0.243465D 05 RHIN= 0.654738D 04 RHAX= 0.421457D 05 ENERGY 3 -0.818600D 01 
PERIOD 3 0 .378065D 05 HHAG= 0.672142D 05 H VECTOR 0.2352730 04 0.240013D 04 -0.671301 D 05 
BNAG= 0.731076D 00 E VECTOR 0.503439D 00 -0.530113D 00 -0.130912D-02 BMAG= 0.421454D 05 
BURN ARC 

LEG 3 4 STATE AT END -0.289233D 05 0.306526D 05 0.197156D 03 0.165880D 01 0.152538D 01 0.212397D 01 

COSTATE AT END 0.397374D 00 0.348802D 00 0.163531D 01 0.3518320-04 -0 . 363556D- 04 0.3419150-06 

HASS AT END OF LEG 3 0.413055D 07 NUHBER OF INTEGRATION STEPS IN BURN ARC 3 2 

TARGET VEHICLE FINAL POSITION AND VELOCITY 

R -0.289575D 05 0.3065220 05 0.194053D 03 V 0.161266D 01 0. 150997D 01 0.213822D 01 

RESULTANT ORBIT SPECIFICATIONS 

SBHIBAJOR AXIS 3 0.427401D 05 SHIN- 0.420893D 05 RMAX= 0.433908D 05 ENERGY 3 -0.466309D 01 
PERIOD 3 0.879353D 05 HHAG= 0.130508D 06 H VECTOR 0.648045D 05 0.617591D 05 -0.9U9653D 05 
EBAG 3 0. 152255D-01 E VECTOR -0.622105D-02 0.132014D-01 0.434007D-02 RHAG 3 0.421447D 05 

TOTAL BURN TINE 3 0.380359D 03 ARC TINES 0.121185D 04 0.255374D 03 0.187264D 05 0.124985D 03 

DC -0.3429720 02 -0.407428D 00 -O'. 310306D 01 -0.461416D-01 -0. 154 1540-01 0. 142555D-01 

-0. 157282D-05 -0.905817D-02 0.326442D-05 0,0 0.0 9.0 

DETSRHINANT OF E= 0.596844D 10 DIAGONAL OF E -0. 218607D 05 -0.243782D 04 0.163160D 03 -P. 1660051) 05 
-0. 1 1 3989 D 05 -0.253055D 04 -0.326874D 00 -0. 812435D-04 -0 . 9 19880D-0 2 0. 586795D-05 
DU= 0.6890150-02 DUDT= 0.113840D 00 CK= 0.100000D 01 EVT= 0.100000D-07 
CK=NIN OF 0.290270 02 0.53545D 02 0.10000D 01 0.32418D 04 0.436360 04 0.339B2D 04 

CHANGE REQUESTED IN INITIAL COSTATE, SWITCHING TINES AND FINAL TIME 
-0.36320 3D-0 2 -0. 177758D-02 -0 .557877D-0 2 0.359956D-05 -0 .292653 D-05 -0. 3 14151D-05 

-0.294391D 00 -0.255003D 00 0.189075D 01 0.187236D 01 

END OF ITERATION NUHBER 3 

NEW QO 0.4150260 00 -0.909741D 00 0.111645D-01 -0. 467025D-03 -0.110953D-02 -0 . 398306D-0 3 

NEW SWITCH TINES 0.214556D 04 0.240097D 04 0.211295D 05 0.212545D 05 



ITERATION NUMBER 


COAST ARC 

LEG= 1 STATE AT END 0.524198D 04 -0.390539D 04 -0.665058D 01 -0.467995D 01 -0.625804D 01 - 0. 3 1 1O61D-0 1 

COSTATE AT END -0.105469D 01 -0. 131068D 01 -0.332243D 00 -0. 130879D-02 0.886090D-03 -0. 644332D-0 4 

PSY= 0. 185056D 00 ALPHA= -0.608892D 02 CALCULATED COAST TIME= n.121156D 04 

SBSIMAJOR AXIS= 0.6S4634D 04 RMIN= 0.653122D 04 RMAX- 0.656147D 04 ENERG Y~ -0.304446D 02 
P0aiOD= 0.527119D 04 HHAG= 0.510820D OS H VECTOR 0.798618D 02 0.194182D 03 -0.S10815D 05 

EBAG= 0.231018D-02 E VECTOR 0.831756D-04 -0. 230866D-02 -0 .864613D-05 RHAG= 0.653685D 04 
BORN ARC 

LEG= 2 STATE AT END 0.362975D 04 -0.550059D 04 -0.697346D 02 -0.808575D 01 -0.624167D 01 -0. 507672D 00 

COSTATE AT END -0.133687D 01 -0. 102644D 01 -0.333130D 00 -0 . 863374D-0 3 0. 129603D-02 0.569163D-04 

BASS AT END OF LEG= 0.692738D 07 NUBBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST ARC 

LSG= 3 STATE AT END -0.291243D 05 0.304936D 05 0.704947D 02 0. 114296D 01 0.110832D 01 0.798648D-01 

COSTATE AT BHD 0.358181D 00 0.339522D 00 0.164611D 01 0.334510D-04 -0. 347154D-04 0.458234D-06 

PSY? 0.757600D 00 ALPH A= -0. 163646D 02 CALCULATED COAST TIHB= 0.187286D 05 

SBMIMAJOE AXIS= 0.243576D 05 R BIN= 0.654736D 04 RMAX= 0.421678D 05 ENERGY= -0.818228D 01 
PERIOD= 0.378323D 05 HHAG= 0.672165D 05 H VECTOR 0.235724D 04 0.240658D 04 -0. 671321D 05 

EM AG= 0.731198D 00 E VECTOR 0.503537D 00 -0.530188D 00 -0. 132549D-02 RB AG= 0.421674D 05 
BURN ARC 

LEG= 4 STATE AT END -0.289545D 05 0.306550D 05 0.198022D 03 0.161287D 01 0.1S0952D 01 0.213770D 01 

COSTATE AT END 0.362347D 00 0.335170D 00 0. 164610D 01 0.332106D-04 -0.349383D-04 -0. 635341D-06 

BASS AT END OF LEG= 0.413008D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 2 

TARGET VEHICLE FINAL POSITION AND VELOCITY 

R -0.239545D 05 0.306S50D 05 0.198056D 03 V 0.161295D 01 0.150966D 01 0. 21382-2D 01 

RESULTANT ORBIT SPECIFICATIONS 

SEMIMAJOB AXIS= 0.421548D 05 RHIN= 0.421416D 05 RHAX= 0.421679D 05 ENERGY= -0.472784D 01 
PERIOD= 0.861 351 D 05 HSAG= 0. 129626D 06 H VECTOR 0.652323D 05 0.622156D 05 -0.931500D 05 

EM AG= 0. 312659D-03 E VECTOR 0.223096D-03 -0. 218821D-03 0.100806D-04 RHAG= 0.421679D 05 

TOTAL BURN TIME= 0.380380D 03 ARC TIBBS 0.121156D 04 0.255413D 03 0.187286D 05 0.124965D 03 

DC 0. 1449 16D-01 0.445962D-01 0.346548D-01 0.830917D-04 0.140930D-0 3 0 . 520 9 7 50- 03 

- 0,377 93 3D— 08 -0 . 432389D-03 0.520904D-07 0.0 0.0 0.0 

DETERMINANT OF E= 0.598783D 10 DIAGONAL OF E -0.217546D 05 -0.245682D 04 0.142653D 03 -0.1918O7D 05 
-0.113223D 05 -0.251065D 04 -0.325230D 00 -0. 90 1 942D-04 -0 ,82885 1 D-02 0.592425D-05 
DU= 0. 1 37203D-03 DUDT= 0. 174784D-02 CK= 0.100000D 01 EVT= 0.100000D-07 
CK=HIN OF 0.14577D 04 0. 34878D 04 0.10000D 01 0.98475D 05 0. 14340D 07 0.24517D 04 

CHANGE REQUESTED IN INITIAL COSTATE, SNITCHING TIMES AND FINAL TIME 
-0.354769D-04 -0. 1779640-04 -0. 131337D-03 0.308003D-07 -0.333872D-07 0.730U10D-07 

-0 .2286 85 D-02 -0 . 358370D-02 -0 . 1 01 139D-01 0.153723D-01 

END OF ITERATION NUMBER 4 

NEH QO 0. 4 1499 ID 00 -0.909759D 00 0.110331D-01 - 0. 466994D-03 -0. 1109560-02 -0 . 398233D-03 

NBH SWITCH TIMES * 0.214555D 04 0.240096D 04 0.211295D 05 0.212545D 05 



ITERATION NUMBER 5 


COAST ARC 

LEG= 1 STATE AT END 0.524 199D 04 -0.390538D 04 -0. 6650510 01 -0. 467993D 01 -0.625&O5D 01 -0 311^610-01 

COSTATE AT END -0.105476D 01 -0. 1310740 01 -Q.332199D 00 - 0 . 1 30885D- 02 0. 8861030-03 -0 . 6427000-04 

PSY= 0. 185056D 00 ALPHA= -0.608892D 02 CALCULATED COAST TIUB= 0.121156D 04 

SEMIMAJOR AXIS= 0.654634D 04 RMIN= 0.653122D 04 RMAX= 0.656147D 04 ENEBGY= -0.304446D 02 
PERIOD= 0.5271 19D 04 HMAG= 0.510820D 05 H VECTOR 0.798618D 02 0.194182D 03 -0.510815D 05 

E«AG= 0 . 23 10 18D-02 E VECTOR 0.831756D-04 -0. 230866D-02 -0.864613D-05 RMAG= 0.653685D 04 
BURN ARC 

LBG= 2 STATE AT END 0.362977D 04 -0.5500570 04 -0.697216D 02 -0. 808574D 01 -0.6241690 01 -0.507552D 30 

- COSTATE AT END -0.133695D 01 -0.1026490 01 -0.333047D 00 -0 . 863428D-03 0.129607D-02 0.5705580-04 

MASS AT END OF LEG— 0.692741D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC= 3 

COAST AfiC 

LEG= 3 STATE AT END -0.291243D 05 0.304936D 05 0.7051060 02 0.114296D 01 0 110833D 01 0.7984600-01 

COSTATE AT END 0.3580790 00 0.3394770 00 0.164S77D 01 0.334458D-04 -0.3471180-04 0.4344090-06 

PSY= 0. 757600D 00 ALPH A= -0.163646D 02 CALCULATED COAST TIHE= 0.187286D 05 

SEHIMAJOR AXIS= 0.243576D 05 RMIH= 0.654736D 04 RMAX= 0.421679D 05 ENERGY^ -0.818228D 01 
PBBIOD= 0. 378323D 05 H«AG= 0.672165D 05 H VECTOR 0.235665D 04 0.240605D 04 -0.671321D 05 
EM AG= 0.731198D 00 E VECTOR 0.503537D 00 -0.530189D 00 -0 . 132577D-02 RBAG= 0.421674D 05 
BORN ARC 

LEG= 4 STATE AT END -0.289545D 05 0.306551D Q5 0.198089D 03 0.161295D 01 0.150966D 01 0 213822D 01 

COSTATE AT END 0.362244D 00 0.335124D 00 0.164576D 01 0.332053D-04 -0. 349747D- 04 -0. 659162D-06 

MASS AT END OF LEG= 0.412954D 07 NUMBER OF INTEGRATION STEPS IB BURN ARC= 2 

TARGET VEHICLE FINAL POSITION AND VELOCITY 

R -0 . 239545D 05 0.306551D 05 0.198089D 03 V 0.161295D 01 0.150966D01 0.213822D 01 

RESOLTANT ORBIT SPECIFICATIONS 

SEHIMAJOR AXIS= 0.421678D 05 RMIN= 0.421677D 05 RMAX= 0.421679D 05 ENERGY= -0.472637D 01 
PERIOD= 0.861752D 05 HMAG= 0.129646D 06 H VECTOR 0.652483D 05 0.622307D 05 -0.931S67D 05 
EfUG= 0.271196D-05 E VECTOR 0.171975D-05 -0 . 208828D-0 5 - 0 . 1 90480D-06 R8AG= 0.421679D 05 
TOTAL BURN TISB= 0.380404D 03 ARC TIMES 0.121155D 04 0.255412D 03 0.187286D 05 0.124992D 03 

DC — 0 . 384675D— 03 0.317517D-04 -0.196088D-04 0.131628D-06 0 . 1 73770 D-07 -0. 174528D-07 

0.14842SD-10 0. 460771D-06 -0.24 0278 D- 10 0.0 0.0 0.0 

DETERMINANT OF E= 0.599584D 10 DIAGONAL OF E -0.217534D 05 -0.245960D 04 0.142443D 03 -0.191B03D 05 
-0.1 133 OSD 05 -0 .251 111D 04 -0.325336D 00 -0 . 903365D-04 -0 .827456D-02 0.592797D-0S 
DD= 0.518046D-07 DUDT= 0.182016D-05 CK= 0. 100000D 01 EVT= 0.100000D-07 
CK=MIN OF 0.38607D 07 0.33492D 07 0.10000D 01 0.78409D 08 0.48664D 10 D.99215D 08 

CHANGE REQUESTED IN INITIAL COSTATE, S HITCHING TIMES AND FINAL HME 

0. 3424 14D-07 0.160488D-07 0.354073D-07 -0 . 35564 ID- 10 0.243606D-10 - 0 . 785 175D- 1 0 

0.956533D-06 0.258524D-05 0.660976D-06 0.310746D-07 

. END OF ITERATION NUMBER 5 

NEW QO 0.414991D 00 -0.909759D 00 0.110332D-01 -0. 466994D-03 -0. 1 10956D-02 -r>. 398233P-03 

NEW SWITCH TIMES 0.214555D 04 0.240096D 04 0.211295D 05 0.212545D 05 



ITERATION NUMBER 6 


COAST ARC 

LEG* 1 STATE AT END 0.524199D 04 -0.390538D 04 -0.665051D 01 -0.467993D 01 -0.625805D 01 -0 . 3 1 106 1D-0 1 

COSTATE AT END -0.105476D 01 -0.131074D 01 -0.332199D 00 -0 . 1 30885D-O 2 O.886103D-O3 -0 6U270On-0a 

PSY* 0. 185056D 00 ALPHA* -0.608892D 02 CALCULATED COAST TIME= 0.121155D 04 

SEBIBAJOH AXIS= 0.654634D 04 RBIN= 0.653122D 04 RBA X= 0.656147D 04 ENERGY* -0.304446D 02 
PERIOD* 0.527119D 04 HMAG* 0.510820D 05 H VECTOR 0.798618D 02 0.194182D 03 -0.510815D 05 

EM AG* 0. 2310 18D-02 E VECTOR 0.831756D-04 -0. 230866D-02 -0.864613D-05 RHAG* 0.653685D 04 
BURN ARC 

LEG= 2 STATE AT END 0.362977D 04 -0.550057D 04 -0.697216D 02 -0.808574D 01 -0.62U169D 01 -0.807552D 00 

COSTATE AT END -0.133695D 01 -0. 102649D 01 -0.333048D 00 -0 . 863 428D-0 3 0. 129607D-02 0.570588D-04 

MASS AT END OP LEG= 0.692741D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC* 3 

COAST ARC 

LEG* 3 STATE AT END -0.291243D 05 0.304936D 05 0.705106D 02 0.114296D 01 0.110833D 01 0.798460D-P1 

COSTATE AT END 0»358079D 00 0.339477D 00 0.164577D 01 0.334458D-04 -0 . 347 1 18D-04 0.4344160-06 

PST* 0. 757600D 00 ALPHA* -0.163646D 02 CALCULATED COAST TIHE= 0.187286D 95 

SEMIMAJOR AXIS* 0.243576D 05 RMIN* 0.654736D 04 RMAX* 0.421679D 05 ENERGY* -0.818228D 01 
PERIOD* 0. 378323D 05 HMAG* 0.672165D 05 H VECTOR 0.235665D 04 0.240605D 04 -0.671321D 05 
EMAG* 0-731198D 00 E VECTOR 0.503537D 00 -0.530189D 00 -0 . 132577D-02 RBAG* 0.421674D 05 
BURN ARC 

LEG* 4 STATE AT END -0.289545D 05 0.306551D 05 0.198089D 03 0.161295D 01 0.150966D 01 0.213822D 01 

COSTATE AT END 0.362245D 00 0.335125D 00 0.164576D 01 0.332053D-04 -0.349347D-04 -0 . 6591 55D-06 

MASS, AT END OP LEG* 0.412954D 07 NUMBER OF INTEGRATION STEPS IN BURN ARC* 2 

TARGET VEHICLE PINAL POSITION AND VELOCITY 

R — 0. 289545D 05 0.306551D 05 0.198089D 03 V 0.161295D 01 0.150966D01 D.213822D 01 

RESULTANT ORBIT* SPECIFICATIONS 

SEMIMAJOR AXIS* 0.421678D 05 RMIN* 0.421677D 05 RMAX* 0.421679D 05 ENERGY* -0.472637D 01 
PERIOD* 0. 861752D 05 HMAG* 0. 129646D 06 H VECTOR 0.652483D 05 0.622307D 05 -0.931567D 05 
EMAG* 0.266009D-05 E VECTOR 0.170161D-05 -0 . 20376 4D-0 5 -0 . 1 69356D-06 RMAG* 0.421679D 05 
TOTAL BURN TIME* 0.380404D 03 ARC TIMES 0.121155D 04 0.255412D 03 0.187286D 05 0.124992D 03 



TEST CASE NUMBER 2 FOR DOUBLE PRECISION DECK 
RENDEZVOUS , 44 DEGREE PLANE CHANGE (7/0/69) 
APRIL 28 , 1970 




SUMMARY TABLES 


ITERATION TOTAL BORN TIME 
NOMBER 


LENGTH OF BURN AND COAST ARCS 


1 

2 

3 

0 

5 

6 


0.373 97 79D 03 
0.369 1784D 03 
0. 3803586 D 03 
0. 3803796D 03 
0 . 38040 38 D 03 
0.3804038D 03 


0.11 13687D 04 
0.1201057D 04 
0.121 18 SOD 04 
0.1211555D 04 
0. 121 15S3D 04 
0. 121 1553D 04 


0.2551809D 03 
0.2553634D 03 
0.2553738D 03 
0 . 2554131D 03 
0 . 25541 18D 03 
0.25541 18D 03 


0.1872919D 05 
0. 1876524D 05 
0. 1872644D 05 
0. 1872858D 05 
0. 1Q72857D 05 
0. 1872857D 05 


0 . 1 187970D 03 
0. 1 138150D 03 
0.124 9849D 03 
0. 1 249665D 03 
0.1249920D 03 
0.124 9920D 03 


ITERATION 

NUMBER 


ERROR IN BOUNDARY CONDITIONS- DC (1-8) 


1 

2 

3 

4 

5 


0.251195D 04 
-0 . 30936 2D 03 
-0.342972D 02 
0. 14491 6D-01 
-0.384675D-03 


0.441559D 04 
0. 317683D 03 
-0. 407428D 00 
0.445962D-01 
0.317517D-04 


-0. 172545D 03 
0 .1 95678D 02 
-0.310306D 01 
0.346548D-01 
-0.196088D-04 


0. 538835D 00 
0. 135471D 00 
-0.461416D-01 
0.830917D-04 
0. 1 31 628D— 06 


0.746112D 00 
0.159980D 00 
-0. 154154D-01 
0. 1409 30D-0 3 
0 . 173770D-07 


0.76281BD-01 
0 . 188098D 00 
0.142555D-01 
0.520975D-03 
-0,1745 28 D-07 


-0. 291278 D-03 
0. 176786D-04 
-0. 15728 2D-05 
-0.377933D-0R 
0. 148425D-1 0 


-0. 1265 1 2D 00 
-0. 857 1 54D-01 
-0.905817D-02 
-0. 432389D-03 
0. 460771 D-06 


ITERATION DC (9-12) 

NUMBER 


DU 


DUDT CK 


0.191551 1 D-04 0.0 

0*0 

0.0 

-0. 2244607D-04 0.0 

0.0 

0.0 

0 . 32644 18D- 05 0.0 

0.0 

0.0 

0.5209039D-07 0.0 

0.0 

0.0 

-0. 2402782D-10 0.0 

0.0 

0.0 


0.1012535D 00 
0.5553624D-01 
0 . 609 01 47D-02 
0.1372028D-03 
0.5180462D-07 


0. 1 657 93 8 D 01 
0. 1 029662D 01 
0. 1138397D 00 
0. 1747843D-02 
0. 1820162D-05 


0. 1000000D 01 
0.100O000D 01 

o. loonoooD oi 
o. lononooo oi 

0.10000003 01 


ITERATION 

NUMBER 


NEW QO GENERATED BY PRESENT ITERATION 


0 

1 

2 

3 

4 

5 


0 . 460 17880D 00 
0. 40157417D OO 
0, 4 1866804D 03 
0.4 1 5026 16D 03 
0.41499068D 00 
0. 41 499071D 00 


-0. 88685450D 00 
-0. 9151 3656D 00 
-0. 90798498D 00 
-0. 90974098D 00 
-0. 90975876D 00 
-0.90975875D 00 


0 .4 1527300D- 0 1 -0 . 43940000D-03 -0 . 1 0504000D-02 -0 . 408 100 00 D-0 3 
-0.3 5542 380D-01 -0.45371994D-03 -0. 1 1 160821D-02 -0 . 44260719D-03 
0 . 16743494D-01 -0. 47063562D-03 -0 . 1 1066283D-02 -0 . 39517370D-03 
0.1 1164458D-0 1 -0 . 46702498D-03 -0 . 1 1 095 285D- 0 2 -0 . 39 830576D- 03 
0.11033122D-01 -0.46699417 D-03 -0 . 1 1095618D-02 -0.39823272D-03 
0.1 1033157D-0 1 -0.466 9942 ID- 03 -0 . 1 10956 1BD-0 2 -0. 39823280D-03 



